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Estimation of Diffusion Coefficient 
of Lanthanum Ions from One-Dimensional 
Liesegang Formation 


K V KURIENT, V K VAIDYAN & M A ITTYACHEN* 
Department of Physics, University of Kerala, Kariavattom, Trivandrum 695 581 


Received 5 July 1982 


Experimental conditions for obtaining Liesegang rings in the growth of lanthanum molybdate crystals in silica gel by single 
diffusion technique have been discussed. The time law, spacing law and the law relating diffusion depth and width are verified. 
Making use of Fick’s law of diffusion, a method has been devised to estimate the diffusion coefficient of lanthanum ions in silica 
gel and the various parameters under which the coefficient varies are investigated. 


1 Introduction 

Diffusion is the spontaneous equalization of 
concentration of molecules, ions or colloidal particles 
in a system as a result of their thermal chaotic motion. 
The role of diffusion in crystal growth has been 
reviewed by several authors'**. The work of Kirov? 
could elucidate the peculiarities of the diffusion and 
crystallization which help to control the growth of 
crystals by diffusion methods. Making use of a simple 
mathematical model, he could deal with the diffusion 
controlled crystal growth in a qualitative way and 
could calculate diffusion coefficient from the first and 
subsequent precipitates produced in the diffusion 
media due to counter-diffusion of two reacting 
substances. Malton and Packter* have shown that 
solprotection (by adsorption) and diffusion are 
important factors in rhythmic precipitation. Packter 
and Roy” have calculated the diffusion coefficient of 
alkaline earth metal ions in melt and found that the 
values are inversely proportional to their ionic radii. 
Diffusional characteristics of a wide number of 
inorganic ions in agar gel have been studied by Lee and 
Meeks®. Making use of the theory of one-dimensional 
Liesegang phenomena, Shinohara’ could estimate 
diffusion coefficient of the diffusiong ions. 

) The authors have developed a method to estimate 
the effective diffusion coefficient of outer electrolyte in 
gel media, making use of the simple diffusion theory 

and nucleation conditions for the formation of one- 

_ dimensional Liesegang phenomena®. The present 

_ paper reports the Liesegang ring formation of 

hanum molybdate (La,(MoO,),;-XH,O) crystals 
® technique and a method of 
ion coefficient of lanthanum ions. 


2 Experimental Details 

Silica gels were prepared in test tubes using sodium 
silicate solution of specific gravity 1.03 with dilute 
HNO, . A standard solution (1.0 S) was prepared using 
HNO, (NH4),NO,; and (NH,),MoQ, as follows. 
Two grammes each of (NH4),MoO, and (NH,),NO, 
were dissolved in distilled water and to this solution 6 
ml of concentrated NHO, was slowly added and the 
whole system was kept cooled. Finally this solution 
was diluted with distilled water to make the total 
volume 100 ml. 

The (NH4),MoO, impregnated in the gel acted as 
the inner electrolyte. The acid content in the standard 
solution was utilized to adjust the pH of the gel 
medium. A set of test tubes with gelling mixtures of 
different pH values between 6 and 7.5 were prepared. 
After the gels were properly set, the outer electrolyte 
(lanthanum nitrate) solution of particular con- 
centration was poured into the test tube. As soon as the 
outer electrolyte was poured a thick precipitate band 
was formed below the gel interface, below which 
Liesegang rings were formed. The distance between the 
successive bands was found to increase with the 
distance of the band from the interface and the 
concentration. The following rules: 


X,/t}'? =a constant (Ref. 10) mops (if) 
X,+1/X, =a constant (Ref. 11) ox 042) 
and 

X/AX =a constant (Ref. 12) cael p 


were found to be obeyed, where X,, is the distance of the 
nth ring from the interface, ¢, the corresponding time of 
formation of the ring, and AX is the distance between 
successive rings. Table 1 gives the parameters of 
Liesegang rings of a typical lanthanum molybdate 
system (Fig. 1). 
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Fig. |—A typical Liesegang ring system of 
La,(MoO,)XH,0O in silica gel medium 


Table 1—Parameters of Liesegang Rings of Lanthanum 
Molybdate Crystals 
[Experimental details: Silica gel of density 1.03 gem *: pH of gel: 


6.5; Outer electrolyte : 0.5 M La(NO;),.6H,O; Age of gel: one day: 
Inner electrolyte : 0.25 S(NH4).MoO,: Working temperature range: 


38-32 C] 

; Xn Eas 
Ring X er: 9 | Milt KPAX A" — 
No. x, tn 
(nr) © br cm cm?@hr! 

J 5.185 68 0.153 0.395 33.889 1.050 1.050 

+ S717? 8 0168 0.403 34.030 1.050 1.060 

7 6271 99 0.184 0.397 34.082 1.050 1.040 

9 6913 117 0.203 0.408 34.054 1.050 1.054 


ek 


3 Estimation of Diffusion Coefficient 
In order to determine the effective diffuison 
coefficient, it is assumed that C4, the concentration of 
the outer electrolyte at distance X Irom the gel 
oe f seconds after the diffusion Starts, is given 
y: 


Ca = Coa exp(— X7/4D 41) 


... (4) 
the outer electrolyte 
sion coefficient of A 
|. Taking a large 
assume that Co, is 


ce-X from the gel 
Position, so that 
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Cz = Co,exp(— X7/4D at) 6) 
Hence the effective diffusion coefficient of A ions, which 
is the result of the movement of both A and B ions in 


the medium, is 


D,= x?*/t ...(6) 
4 In (Co4/Cp) 
Thus by using Eq. (6), D4 is calculated. 

In the actual experiment, three arbitrary values of 
inner electrolyte, viz. 0.25, 0.5 and 1.0 S, corresponding 
to Cz values 1.327, 2.666 and 5.310 x 107° ions/cm?, 
respectively, against the outer electrolyte concen- 
tration in the range 0.25-0.7 M were selected. The rate 
of movement of the advancing boundary X and time t 
were noted and the mean values of X?/t (Table 2) were 
calculated and D, values estimated. The uncertainty in 
the determination of D, is of the order of 10%. Within 
this limit, the diffusion coefficient is estimated as (0.73 
+0.07) x 107° cm’s“'. For fixed value of Cg (0.25 S), 
values of D4 were estimated with Co, being varied from 
0.25 to 0.7 M. A small but observable diminution in the 
X?/t and D, values were observed due to C g Variation, 
for a fixed value of Co,. The increase in inner 
electrolyte could cause a small retarding effect on the 
onward moving outer electrolyte. 


Table 2— Estimation of Diffusion Coefficient of Lanthanum 
Ions in Silica Gel 
[Experimental details: Silica gel of density 1.03 gem *; pH of gel 
medium 6.5: Working temperature range 28-32 C] 


Electrolyte Conc. of Conc. of X?/t D, 
A ions B ions (Cg) (cm?.hr') (x 105 
Outer Inner (Cas) ( x 10?°cm~ 3) cm?s_') 
(x 10?°cm~ 3) 
0.5M 0.25S $7.250 1.327 0.400 0.738 
05M 0585 $7.250 2.655 0.325 0.735 
05M 10S 57.250 5.310 0.245 0.717 


Table 3—Effect of Gel Density (gem~ +) and Ageing on X2/t 
(cm*hr') and D, ( x 10°cm?S~') Values 
[Co4 = 57.25 x 10?° ions cm~*; Cg = 1.327 x 10?° ions .cm~} 
Gel Density 
eT 


pH 1.025 1.030 1.035 1.040 
Fe Tse | Per X*t Dy Xt D, 
6.0 0.390 0.730 0.400 0.738 0.402 0.745 0.404 0.750 
6.5 0.392 0.728 0.394 0.727 0396 0.742 0.398 0.745 
7.0 0.389 0.720 0.391 0.722 0.393 0.737 0.395 0.742 
Ageing 
SE 
One day Two days Three days 


6.0 0.410 0.756 0.404 0.745 0.400 0.738 
65 0.400 0.738 0.396 0.730 0390 0.720 
70 0.390 0.720 0.385 0.710 0.380 0.700 


ee 


—— 


ie a. \ fe tl ee 


KURIEN e7 al.: DIFFUSION COFFICIENT OF LANTHANUM IONS 


As a result of mild interaction between gel and 
diffusing ions during the process of diffusion, the 
eflective concentration of the outer electrolyte is 
controlled by the density of gel as derived by Malton 
and Packter*. In the present case, a slight increase in 
the density of the gel medium is found to act as a 
retarding medium for diffusion and retards the 
reaction of precipitation. Table 3 shows the change in 
D, values with gel density variation. The age of gel is 
also found to affect a considerable decrease in the 
estimated values of D, (Table 3). As the density 
increases, the cell wall of the gel thickens” and as the 
age increases, the gel gets shrunk, so that in both cases 
the pore size reduces'*, which in turn decreases the 
ionic diffusion. This explains the observed decrease in 
D, with gel density and ageing. 
eed 

The one-dimensional Liesegang phenomena can be 
utilized to estimate the effective diffusion coefficient of 
outer electrolyte. The diffusivity of the medium 
_ diminishes with increase in gel density and ageing. 
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Effects of Quenching & y-ray Irradiation 
on the Dielectric Properties 
of Calcite Single Crystals 
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ielectri K), loss (tan 0) and conductivity (a) of calcite single crystals irradiated with j-rays before and after 
rae Soe is pone of frequency in the region 107-10’ Hz and in the temperature range 30-400 C with the 
ee optic axis of the crystal. Value of K at 30 S (equal to 7.8 and he 2 
quenching and )-ray irradiation, however, the values of tand, particularly at ow irae aac: - pe x: ‘ 
irradiation and quenching. Increase of K with temperature.at different frequencies is more marked for pray gee pis sample. 
the rate of increase of K at any frequency is maximum for the quenched and }j-ray irradiated sample. A similar aviour is 
exhibited by tan 0 versus temperature graphs for these samples. The plots of log o versus 1/T at different ae eae 
a straight line beyond 300 C; the activation energy values for conduction in this region are found to be 0.9, 0.73 a 0. is 
respectively for the as-cleaved, y-ray irradiated, and quenched and y-ray irradiated samples. An attempt is made to explain these 


observation on the basis of space charge effects. 


1 Introduction 

Natural calcite is available in the form of large single 
crystals which are transparent in the near visible 
region. These crystals show an extraordinary long- 
period after-glow when irradiated with lonizing 
radiations like X-rays. As such, the optical properties 
of both pure and coloured calcite have been extensively 
studied’ ~*. Recently, possible models for colour 
centres responsible for absorption bands in X-ray 
irradiated calcite crystals have been suggested?. 

When X-ray irradiated calcite crystals exhibit 
phosphorescence or thermoluminescence resulting in 
release of charge carriers, itwas observed that the 
dielectric properties of the crystals were considerably 
affected. Such studies, carried out in this laboratory® ’, 


gave useful information regarding the defect processes 
taking place in this solid. 


)-Ray irradiation is expected to produce consider-: 


able changes in the dielectric Properties of calcite 


dielectric 
Properties of natural calcite single crystals before and 


. These roperti . 
with electric field parallel Properties are studied 


known to the authors, such a stud 


reported till now-though the dielectr; . 
; ric pro 
calcite crystals were investigated earlier® emma: 
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2 Experimental Details 

The present work was carried out on cleaved crystals 
of natural calcite of good optical quality. The samples 
were ground and polished and had, in their final form, 
the dimensions of 20 x 20x 1 mm?. The orientation 
was checked by means ofa polarizing microscope. The 
dielectric measurements were taken on a GR 716 type 


capacitance bridge” up to 105 Hz and on a Marconi 


circuit magnification metre type TF 329G in the range 
10°-10’ Hz using the resonance curve principle®. Silver 
paint was applied on both sides of the crystals to serve 
as electrodes. The accuracy of the measurement of the 
dielectric constant K is 2% and the loss (tan 0) is to 
about 5%. Three samples of calcite crystals were 
studied and the agreement was within the limits 
mentioned. Flame photometric and chemical analysis 
of the calcite crystals showed that Sr is present to about 
10 ppm and Mn to 4 ppm. 

The samples were irradiated with )-rays from a® Co 
source of 1000 Ci for 48 hr; the samples were kept 
Practically touching the source. Quenching of the 


crystals to room temperature was done after heating 
them at 800°C for 1 hr. 


3 Results 

The dielectric constant of calcite crystals at 30°C and 
10° Hz with E|JC is measured to be 7.8 and is 
frequency-independent, this value remz — un- 
affected by quenching and )-ray irradiation; however, 
the dielectric loss is considerably increased (Fig. 1) due 
to )-ray irradiation and uenching. . 


VEERAISH er a/.: DIELECTRIC PROPERTIES OF CALCITE SINGLE CRYSTALS 
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Fig. !— Dielectric loss (tan 5) at 30°C as a function of frequency of 
calcite crystals under different conditions [(a), as-cleaved; (b), y-ray 5 0 10 200 300 100 
irradiated; and (c), quenched and )-ray irradiated] Temperature ,T 


Fig. 3—Dielectric constant (K) as a function of temperature at 5 

x 10? and 10° Hz under different conditions for calcite crystals [(a) 

and (a’), as-cleaved; (b) and (b’), y-ray irradiated and (c) and (c’), 
quenched and )-ray irradiated] 
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* Fig. S—Conductivity (¢) of y-ray irradiated calcite crystals as a 
function of 1/T at different frequencies 


constant of a material at low frequencies; electronic, 
tonic, dipolar and space charge. The nature of variation 
_of K with frequency indicates which contributions are 
present. Space charge polarization depends on the 
purity and perfection of the crystals. Its- influence is 
noticeable mainly in the low-frequency region. The 
dipole-orientational effect can be ex hibited even up to 
10'° Hz. The ionic and electronic polarizations always 
exist below 10!3 Hz. 
y-Ray irradiation of a material is known to knock off 
_ tons from their positions in addition to producing 
electrons and holes'°. The displaced ions may collect at 


Table 1—Summary of Data on Dielectric Properties of Calcite Single Crystals 


the surface of the crystal and some of them may also 
occupy interstitial positions. As the carbonate bond in 
calcite is apparently weak, )-ray irradiation of these 
crystals may knock off oxygen ion: ; additionally, some 
of the other ions of calcite lattice may also be removed 
from their positions. So the final effect on ‘y-ray 
irradiated calcite crystals seems to be the generation of 
a large number of charged defects like vacancies and 
interstitials. Quenching and y-ray irradiation of the 
crystals are expected to enhance the concentration of 
these defects. 

Reviewing our data, the considerable increase in 
dielectric loss at room temperature (=~ 30°C), 
particularly at low frequencies, exhibited by calcite 
crystals irradiated with y-rays, before and after 
quenching, may be ascribed to space charge 
polarization due to the defects. However, the 
maximum value of tand measured at 10? Hz in 
quenched and )-ray irradiated calcite is only 0.01 
which is still small enough not to affect the dielectric 
constant at that frequency. As such, the quenched and 
)-ray irradiated samples did not show any increase in 
their K values at room temperature. 

Temperature has a complicated influence on the 
dielectric constant. Increasing the temperature 
generally decreases the electronic polarization slightly. 
The increase in ionic polarization due to thermal 
expansion is also not very large. Even assuming the 
Presence of a small number of dipoles and their 
contribution to dielectric constant, we know from 
Debye’s theory that K is inversely proportional to 
temperature. As such, it is expected that the dielectric 
constant of calcite should not change considerably 
with increase in temperature. But we find large increase 
in K values with temperature in calcite crystals )-ray _ 
irradiated before and after quenching. The crystal | 
defects only can increase and also probably become 
more mobile with temperature making the space e 


As-cleaved j-ray irradiated Quenched and y-ray 
irradiated 7, 
Para- F 7 i OE Os ee “eg ae 
meter <7 Sn WP Bn ee 10? Hz 10° Hz 10? z 10° Hz oe 
hla cada od Ee ee Ee = 300 150 300 30 150 300 
K ‘ rT 
ip esa ee a eC 78 82 98 78 45 — 78 g8 
tand 20x 25x 10x Lix tix 14x 30x 80x 45x 18 1s ; vs eas 
10°* 10°* 10-3 49-4 40-4 10-* 1-7 40-! as x LSx 10x 95x Six 25x 40x 
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charge polarization dominant: hence K increases with 
temperature''''?. The changes in K with temperature 
in these samples are smaller at higher frequency as this 
type of polarization decreases considerably at such 
frequencies. Our results show that the rate of increase 
of K with temperature (at any frequency) is largest for 
the quenched and y-ray irradiated sample, indicating a 
correspondingly largest concentration of defects in 
that sample. This conclusion seems to be supported by 
the fact that the activation energy for conduction in the 


intrinsic region in these samples is appreciably 
decreased '? (Table 1). 
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Analytical expressions for the thermal noise in single injection solid state diodes with traps lying above the Fermi level have 
been evaluated by using regional approximation and salami methods. The mobility of the injected current carriers is 
proportional to the carrier density. The thermal noise is shown to be highly suppressed at high injection level of current. 


I Introduction 

It is well known that the main origin of noise in 
space-charge-limited current flow through insulators is 
thermal noise’ ° which is generated due to the 
scattering of injected current carriers with the lattice 
imperfections present in the device. The presence of 
traps in the insulator drastically affects the current- 
voltage characteristic. Therefore, the traps are 
influencing the thermal noise generated in the device. 
The noise calculations may be easily obtained by the 
Salami method given by Van der Ziel':*. In this 
method, the insulator is divided into small slabs which 
are having noise source present within them. 

In a low mobility insulator, the mobility of the 
current Carriers is proportional to carrier density. The 
results for space-charge-limited single injection current 
in low mobility insulators are obtained with the help of 
regional approximation method®~ 8. In this paper, the 
thermal noise is studied to show the effects of space 
charge, traps lying above the Fermi level and carrier 
density dependent mobility. 


2 Regional Approximation Method 
: Let us consider a current injection in a low mobility 
insulator containing a significant density of a single set 
of traps lying above the Fermi level. The general 
equations governing the single injection current flow 
and Poisson’s law in non-constant mobility regime are 
given by: 

J =epnE .. (1) 
&\dE 
Gee ey Pe 


where x is the mobility of the current carrier, n and No 
are the concentrations of free electrons and their 
thermal equilibrium value and n, and No are the 


concentrations of trapped Carriers, and their thermal 
equilibrium value respectively, 
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Wintle® has shown that in a low mobility insulator 
containing trapping sites, there is no delay to the 
moving carriers due to the trapping effect, because the 
trapping states are permanently occupied by electrons. 
It gives rise to the mobility (4) of the current carriers 
proportional to the carrier density: 
ft = hn(x) Bere Oo 
where h is the proportionality constant. These 
equations are subjected to a boundary condition for 
ohmic contact 
E(0)=0 ...(4) 

At an appropriate injection level of current, the 
insulator is divided into four separate regions with the 
help of the regional approximation method>7°8, These 
regions are well separated by three imaginary 
transition planes x,, x, and x3 which are shifted 
towards the anode with the increase of injection level of 
current. The three regions I-III are space-charge- 
limited and the region IV is ohmic. Agrawal et al.® have 
given the detailed description of the application of the 
regional approximation method to the present 
problem. The results for transition planes are given by® 


Sane (5) 
1 3e7hB ns = 


ee (6) 
2 hBCns af 


20eJ 


“3 3ethnd ee fe 


where B = N,/no, C = N/gno and 6 = C/B -«. (8) 
here N, is the total number of traps, g the statistical 
weight of traps and N has already been defined 
elsewhere’. It is obvious from Eqs (5)-(7) that the 
transition planes are proportional to the current 


density. Therefore, the planes are reaching towards the 
anode with increase of current. 
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The complete current- -voltage characteristic is 
divided into four separate regimes by three critical 
Currents defined as 


M3 V= Lh, X22) =L, X,VSor3)=L ...(9) 
3 Thermal Noise Calculations 

The scattering of the current carriers with the lattice 
imperfections of the insulator yields the current or 
voltage fluctuations which is the main source of noise 
in the space-charge-limited solid state diode. A very 
simple method of noise calculations in such a device is 
given by Van der Ziel*. In this method, the insulator is 
divided into slabs of small thickness by the planes 
perpendicular to the current flow. The noise is 
generated in each slab due to the presence of noise 
sources in them. The total noise voltage across the 
device is obtained by taking the sum of voltage drops 
across each slab. 

According to Van der Ziel*, the noise voltage in a 
frequency interval Af in the open-circuited device may 
be represented as the mean square noise voltage given 


by: 
Av? = 4kT ARAS. AR = R= (10) 


where AR is the small resistance of a small section Ax, S 
is the area of cross-section and Eq. (1) is used to obtain 
the last equality. Summing over all sections of the 
insulator, Eq. (10) yields 


v=4kTRAS, R=ZXAR= Akh) 


als 


, L 
where V= | E(x)dx fi) 
0 
The thermal noise is generated in the complete range 
of current-voltage characteristic in the order of regimes 
as discussed below. 


-3.1 Ohmic Regime (J <J,,. ;), (x3 < L) 


Initially the injection level of current is very low so 
that the injected space-charge does not reach the 
anode. The current flow is carried by the thermally 

d free-carriers. Therefore, Eqs (1) and (12) with 
| me Ns the current-voltage characteristic in ohmic 


e(h3) 


sista ae (R ta) and thermal noise (va) are dacieed from 


...(14) 


At ~ is 


in nature. This regime will terminate at a critical 
current J.,., when the plane x, lies on the anode. 


3.2 Shallow Trap Regime (J... , <J <J,, ;) 

The three regions I-III are present in the insulator. 
The transition planes x, and x, are very close to 
cathode. Therefore, the main contribution to thermal 
noise is Obtained by the shallow trap region II]. From 
the analysis®, the onset values of shallow trap regime 
are given by: 


3e7*hngLl 

Jor = ot ae ..(15) 
_ 9eng L? 

Fo 1 100¢ .. (16) 


where the critical voltage V.,., is corresponding to the 
critical current J,,.,. The cane noise resistance and 
thermal noise are derived from Eqs (11), (15) and (16) 
given by: 


ee: 
Rot Fo = Sgngts --(17) 
Ver, 1 = 4kTR.,., Af ..(18) 


This regime is terminated at the critical current J,,. >. 


3.3 Trap Filled Limit (TFL) Regime 
Ver.2 <J <Jer,3) 


In this regime, the regions I and II are present inside 
the insulator as the transition plane x, leaves the 
insulator at the onset critical values given by® 


Zee 3 
Lee BC ie (18) 
a Beno i? 
Vow, 2 ce Ie ..(19) 


All the traps are gradually filled during the TFL 
regime. The critical noise resistance and thermal noise 
corresponding to Eqs (18) and (19) are evaluated as: 


on Ver,2 z. L 
Rea = JOS 2CtehndS “ 
Yir,g = AKT Ro, 2 Af ..(21) 


where Eq. (11) is used. 
All the traps are completely filled at the terminating 
current and voltage® 


3B? Le*hn3 


Jer,3 = py; ..(22) 
_ 9eny BL? 
A ee gem ..(23) 


The corresponding critical noise resistance and 
thermal noise is derived from Eqs (11), (22) and (23) 


L.A Cae...” 
Res = SS ~ SBF ehngs AB 
Wa =4kT Rey AS 25) 
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4 Discussion . ee: 
Complete variation of noise resistance [v?/4. ‘TALUS 
with applied voltage is shown in Fig. | on a log-log 
———---— scale by using Eqs (13)-(28). The continuous curve is 
started from a miximum value Rg in the ohmic regime. 
SMALL OW. TRAP The noise resistance in the shallow trap regime (V,,, 
FREE <V<V.,,5) is varied from R,,,, to Re,» which < 


lower than the ohmic resistance. There is a steep fall 
from R,,,, to R.,,3 during trap-filled-limit(TFL)regime __ 
(Vi..2<V <V.,,3). In the last space-charge-lim 


c 


TFL 
ERIE Eqs (13)-(28) are used to draw Fig. 1: 
B= 10°, C=500,h=107'29 meV ‘5° 4 
L=107°m, : . = 
no = 102-3, 2 = 10-9 S— 10 en 
It may be pointed out from Fig. 1 and Eq. (2 
the noise resistance during shallow tra 
regimes are parallel to each other. T 
thermal noise in shallow trap regime is 
as the square of the applied voltage. 
oe 3 : range of current-voltage characteristic 
g. 1 —Log-log plot of noise resistance (:?/4k TAf) versus applied inferences may be drawn: 4 , a 
be be (a) The space-charge highly ppr S : he therm 
rap-Free Regime (J >J.,, 5) _ noise from an initial high vz = ns 


[v? 76 kTAt ] ohm 
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Tetramethylammonium tetrachlorozincate (TMAZnCl,) shows six successive phase transitions in the temperature range 

80-300 K. Polarized Raman spectra of oriented single crystals of TMAZnCl, and the infrared spectra by the KBr pellet 

technique have been recorded from 298 K to 130K, covering all the known phase transitions. An increase in the frequency of 

the Zn-Cl stretching vibration is found in Raman and infrared spectra as the temperature is lowered. Three bands are found in 

the Raman spectrum in the region where CH, frequencies usually occur and, in a similar way, three bands appear in the 

‘ infrared spectrum at low temperature due to the deformation vibrations of the CH; group. The triggering mechanism of the 
phase transition at low temperature can be due to a coupling of the distortion of (ZnCl4)?~_ groups and the CH deformation, » 

where the motion of relatively free methyl rotations is slowed down, favouring one equilibrium position below 161 K. The 

phase transition appears to be an order-disorder one with methyl rotational disorder in the high temperature paraelectric 


phase. 


1 Introduction spectra structure correlation, the possible triggering 
Tetramethylammonium tetrachlorozincate abbre- prechehism, 

viated as TMAZnCl, belongs to the family of crystals : : 

of ‘the type A,BX (A=Rb, Cs, NHq..., B=Zn, Co, 2 Experimental Details eh: 

§ V J i, &=Cl, Br, I...) which’ share an incom- Crystals of [N(CH 3)4],ZnCl4 were grown byaslow | 


eer ne 


3 _ mer isurate phase in their phase transition sequence. evaporation of an aqueous solution containing 
| te nCl, cans pee six phase transitions in the stoichiometric ratios of TMACI and ZnCl,. The 
_ temperature range 80-300 K. They are shown in Table _ crystals were transparent and colourless. Cosa 
See Pach X-ray crystal structure’, in the form of a Jiadon 488 on were ea yy an 
ress and neutron bei studies® Act ion le poe ee ) a? = 
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superstructure which is characterized by a wave vector 
q=(1 —0d)a*/3 where 6 depends on the temperature, 0 
= —0.34 at 7;,=293K and 6=0 at 7,=279.2K. The 
crystal structure is similar to CsZnCl, where Cs~ ions 
are replaced by N(CH;)4 ions'’. There are two kinds 
of polyatomic radicals in the unit cell, one being the 
cation N(CH;); .(8 per unit cell) and the other 
(ZnC1,)* (4 per unit cell). According to group theory, 
TMAZnCl, has 468 modes separated into 396 internal 
modes and 72 external modes (36 rotatory lattice 
modes and 36 translatory lattice modes of which 3 
form the acoustic branches). The symmetry 
correlations for [N(CH;),],ZnCl, in the paraelectric 
and ferroelectric phases are shown in Tables 2 and 3 
respectively. The ions occupy only general positions in 
the other phases below 276 K. 


3 Results 
3.1 Raman and IR Spectra at Room Temperature (298 K) 


Polarized Raman spectra of the single crystal of 
TMAZnCl, have been taken for six orientations in the 
range 50-3400cm~' at 298K and the infrared 
spectrum of this compound was taken in the range 200- 


Table 2—Symmetry Correlation for Paraelectric 
[N(CH,)4], ZnCl, at Room Temperature (293 K) 


Free ions Sites Crystal Activity 
NicH,It 

eisciihy Cs Ooh 
Zncl, 


Raman 

Raman 

Raman 

Raman 

Inactive 

Bry IR 

IR 

| B;, IR 

7) 2 [eee nee 

Table 3—Symmetry Correlation 
[N(CH),], ZnCl, 


Free ion Sites Crystal Activity 
‘4 bd Coy 
A 


for Ferroelectric 


. Ay Raman & IR 
A Ay Raman 
F 
‘ 8 
; Raman & IR 


F, 


A eranme ee ct 
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4000 cm ~'!, using the KBr pellet technique at 298 K. 
The observed Raman and IR frequency shifts at 298 K 
and at low temperature and the assignments are shown 
in Table 4. The typical Raman spectra are shown in 
Fig. |. The polarized Raman spectra of TMAZnCl, 
contains the lattice modes and ZnClz~ internal 
vibrations (50-350cm~'), the ammonium skeletal 
vibrations (350-1050 cm ~') of the TMA®* ions, the 
methyl bending and rocking vibrations (1050- 
1500cm~') and the methyl stretching vibrations 
(2500-3200 cm ~'). In the lattice region, the Raman 
spectrum is generally weak in intensity excepting a 
sharp line at 277cm~' which is interpreted due to 
ZnClj~ internal vibrations. The C—N stretching 
frequencies at 757 cm ~' and 954 cm ~™ are strong in 
intensity and are highly polarized. Both methyl 
symmetric deformation bands are weak but both 
methyl asymmetric deformation bands are strong. 


3.2 Low-Temperature IR and Raman Spectra 


The Raman and IR spectra were taken at low 
temperatures covering all phase transitions listed in 
Table |. The results are presented in Table 4. When the 
crystal is cooled to the ferroelectric phase, the ZnCl2 ~ 
Stretching frequency (1',) is increased from 277 to 
279 cm ~'. In the infrared spectrum, this correspond- 
ing line is shifted from 270-277cm~'. The v> 
symmetric deformation vibration of CH3, which gives 
a weak band at 298K, becomes sharper in the 
ferroelectric phase in the Raman spectra. The y, 
symmetric stretching of CH ; band 2987 cm ~' at 298K 
is shifted to 2988 cm! which is reflected in IR spectra 
also. When the crystal is cooled to 133 K, covering the 
transition from phase V to VI, the v, stretching band of 
ZnClz increases from 277 cm! to 282cm~'! in 
Raman spectrum and from 270 to 285cm~— in IR- 
spectrum. The v, and v4 stretching vibrations of C — N 
increase from 757 to 760 cm ~! and 954 to 956 cm~! 
respectively. The v4 asymmetric deformation of CH; i 
increases in intensity while in the region of v, : 
symmetric deformation of CH 3, at 1420cm~', three 
bands appear as 1408, 1414 and 1428cm' 
respectively in the Raman spectra (Fig. 2). In the 
infrared spectra, the corresponding band shows 
maxima at 1409, 1419 and 1429 cm ~! respectively (Fig. 7 
3). The v, stretching frequency of CH, at 2987cm-! 
decreases in intensity, but leads to three bandsa 


respectively. Also the rocking mode of 
infrared 1 is split into two 
and 1295 cm~'. The wid 
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Table 4—Raman and Infrared Frequencies in TMAZnCl, in Various Phases 


Assignment 


vy lattice 


ZnClj ~ internal 
vibrations 


Skeletal TMA 
vibrations (C — N) 
¥2 Symmetric stretching CH, (72,) 
va-Rocking (£,) CH; 
ve-Rocking CH, (72,) 
¥2 symmetric deformation CH, (72g) 


¥4 asymmetric deformation (E,) CH; 
Nonfundamental 


Nonfundamental 
v, symmetric stretching (Aig) 


298K 


278K 


263K 173K 

IR Raman IR- Raman 
_- 63 _— 63 
— 70 — 70 
a 126 — 131 
277 279 280 282 
— 37) — 372 
457 457 457 457 
— 757 = 759 
950 954 950 956 
1080 1069 1080 1070 
1175 1175 1175 1175 
1288 1286 1288 1286 
1417 1420 1417 1408 
1414 
1428 
1450 1456 1450 1459 
1482 1480 1482 1480 
— 2930 = 2928 
2928 2959 
2960 2987 2960 2979 
2987 
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INTENSITY (arbitrary units) 


3106 3006 2906 1495 1395 
WAVENUMBER (cm7) 


Fig. 2—Raman Spectra of tetramethylammonium _ tetra- 
chlorozincate at different temperatures 


133K 


Whereas a large 


CH; increases at low temperatures. The single band at 
952 at room temperature in the infrared spectra splits 
up into two bands and this could be due to 
nonequivalent C—N bonds in the low temperature. 
The variation of CH; stretching frequency with 
temperature in the JR spectra is shown in Fig. 4. 


4 Discussion . 
The structure of TMAZnC], is similar to CsZnCl,. 
The ZnCl, units are not perfect tetrahedra but three of 
the ZnCl, bonds have lengths 2.274A, 2.285A and 
2.285 A whereas the fourth, Zn-Cl bond, is larger with 
a value of 2.303A ZnCl, tetrahedra thus appear to be 
distorted in the form of flattened tetrahedra?. In the 
Raman and IR spectra, the observation of three 
frequencies, characteristic of the CH, stretching 
vibrations and three frequencies attributable to the 
deformation vibrations of the CH, groups in the 
infrared and Raman spectra of TMAZnCl, and the 
broadening of the line due to the CH; stretching 
vibrations at 133 K suggest that the motion of the 
deformed CH, groups play an important role in the 
transition from phase IV to VI. The decrease in 
intensity of v, stretching vibrations of the CH roups — 
and the temperature sensitiveness of the | a 
indicated above, enable one to conclude that eee Vit 
strong coupling between the deformation of CH, 
groups and the distortion of the (ZnCl,)? ~ in the low _ 
temperature phase and could be related tc 
ferroelectric behaviour in phase II]. 
_The phase transitions at the low temperatt 
viewed, therefore, as an order-dis 
methyl rotation disorder in the h 
paraelectric phase. The recent c: 
reveal that methyl g 
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Single crystals of p-type SnSe have been grown by simplified Bridgeman method. The crystals ber page hae ss! 

ice px ained are: a=4.3A, b=4.05A and c=11. . The absorption 
diffraction technique and the lattice parameter values obtained See: 

spectrum of the crystal in the range 0.5-1.3 eV has been recorded for two light polarizations (E |la and E||b).A detailed analysis 

of the absorption coefficient curve in the fundamental absorption edge region indicated the existence of indirect allowed 

transitions. The values of indirect energy gaps were 0.94 eV for E||a and 0.891 eV for E||b. It was found that the phonon 

assisted transitions required the participation of a phonon with energy 0.022 eV for E||a and 0.009 eV for E||b at different 


energy thresholds. 


1. Introduction 

The binary IV-VI compounds with Ge, Sn and Pb as 
cations and S, Se and Te as anions form a very 
interesting class of semiconducting crystals. These can 
be divided into three groups according to crystal 
structure. The lead chalcogenides (PbS, PbSe and 
PbTe) crystallize in a cubic (NaCl) structure, are very 
easy to grow and have been thoroughly studied!. GeTe 
and SnTe have a rhombohedral structure at low 
temperatures, converting to cubic at T=0°C for SnTe 
and at T=400°C for GeTe. These compounds have 
also been extensively studied?, though not as fully as 
lead compounds. The remaining four compounds, viz. 
GeS, GeSe, SnS and SnSe have an orthorhombic (space 
group D5?) crystal structure2 and are least studied of 
the three groups. These are in many ways intermediate 
between two-dimensional (layer) structures and three- 
dimensional crystals. SnSe is a fairly typical member of 
this interesting class of isomorphic materials. 

The SnSe crystal is made up by tightly bound double 
layers of Sn and Se atoms stacked along the C-axis as 
shown in Fig. 1. This lattice structure® can be viewed as 
a severely distorted NaCl structure. Each atom has 
three strongly bonded neighbours within its 0 
and three more distant weakly bonded neighb 
of which lies in an adjacent layer. The bonding 
layers is weak, being of Van der Waal’s t 
mca a ao cleaved in the (001) plane. 

fact, Bnt expect SnSe to exhibit the 
typical anisotropy of a layer structure. 


wn layer 
Ours one 
between 
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2. Experimental Details 

Single crystals of SnSe were prepared by the 
Bridgeman technique. Fig. 2 shows the schematic 
diagram of the furnace used. A translucent silica tube of 
5 cm diameter and 50 cm length was wound directly by 
using Kanthal A wire (1.78Q/m) covered with 
insulating beads. The muffle was kept in a central 
rectangular hole constructed from hot face good 
quality insulating bricks. The space between bricks and 
muffle was filled with alumina powder. The complete 
length of the silica tube was divided into two zones (20 
and 30 cm respectively), each zone having its separate 
heating element, power supply and temperature 
controller. The temperature profile of the furnace used 
for growth is also shown in Fig. 2. A double walled 
quartz ampoule of 20cm length and 2.2 cm diameter 
was half filled with stoichiometric amounts of 
powdered Sn and Se, evacuated and sealed properly. 
The tin and selenium used were of 99.999 %, purity 
supplied by Koch Light Ltd, England. The ampoule 
was vibrated for half an hour to ensure mixing of 
elements and held vertically in the upper zone for 2 
days at a temperature above 900°C (the melting point 
of SnSe is 860°C). The ampoule was removed from the 
furnace, vibrated again for half an hour and put again 
in the furnace for 6 more hours in order to obtain 
complete reaction between the crystal components. 
Now, the ampoule was gradually lowered with a speed 
approximately 3 mm/hr using a gear box specially 
fabricated for this Purpose. The melt of SnSe began to 
crystallize from the bottom of the inner ampoule. 
reached the temperature of 650°C, it 
1 day to reduce 
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Electron diffraction pattern of SnSe crystal 


Fig. 3 


and then cooled down to room temperature ata rate of 
20 Chr. The grown crystals were approximately 3 cm 
long and 1.5 cm in diameter. The electron diffraction 
photograph of the crystal (Fig. 3) showed that the 
material was a single crystal. From X-ray diffraction 
studies of the crystal powder, the values of lattice 
constants were found to be a=4.3 A, b=4.05A and c 
=11.62A. Subsequent electrical characterization 
measurements indicated that the specimens were of p- 
type with hole concentration in the range 10'°- 
f0'’ em=° 

Thin crystals (thickness 50-250 0m) for the 
absorption studies were obtained by cleavage, which 
was performed with the aid of cellophane tape. The a 
and b axes are contained in the cleavage plane of the 
specimens and were located using standard Laue 
diffraction technique. No measurements could be 
performed for E||c because the crystal was 
considerably soft and did not permit cutting and 
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polishing. The optical absorption data were obtained 
in the photon energy range 0.5-1.3 eV at room 
temperature (300 K). When the sample thickness was 
appropriate, interference fringes were observed in the 
transmitted beam for both directions of light 
polarization employed (E||a, E||b). The values of the 
optical parameters were determined following the 
method proposed by Hazelwood’ using least squares 
refinement. 


3 Results and Discussion 

Fig. 4 shows the spectral variation of absorption 
coefficient and refractive index for two directions of 
polarization (E||a, E||b). The difference in the 
absorption coefficients for the two directions Is greater 
in the vicinity of the absorption edge, converging to a 
single value as the wavelength increases. A detailed 
analysis® of the power dependence of the absorption 
curve in the fundamental absorption edge region 
revealed the existence of indirect allowed transitions. 
For indirect allowed transitions: 


ahv = A\N (hiv—E,+hw,)? +(1+N,)hv—E, —hw;)*! 
ail) 


where E, is indirect band gap energy, fic; is the energy 
of the phonon and N,,, the number of phonons is given 
by,. 


N,= (exp (hw;/kpT)—1'~' 2) 
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Fig 4 Absorption coefficient and refractive index curves of SnSe 
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2? 


The [x(/v)]' * versus hv curves in SnSe for investigation for 
indirect interband transitions 


Fig. 5 


The first term in Eq. (1) describes processes involving 
phonon absorption, while the second term describes 
phonon emission. Thus a plot of [x(hv)]'/? vs hy is 
expected to resolve into a straight line for each 
independent process. Such a plot for SnSe is presented 
in Fig. 5 from which multiphonon Participation is 
evident. Since each curve shows three knees E,,E, and 
E3, it is safe to assume that two different phonons 
participate in the absorption process with the 
threshold energies spaced symmetrically about the 
energy gap value. As for the forbidden energy gap of 
pure crystal, E, can be determined by extrapolating the 
high energy side of the absorption curve as shown in 
Fig. 5 to x=0. In this way, the analysis of the curve 
showed that E, definitely corresponded to the 
midpoint between E, and E,. The values of E, obtained 
for E||a and E||bare E,,=0.94 eV and E,,=0.891 eV, 
respectively. The phonon energy 4E,— E,)is 0.022 eV 
for E|\|a and 0.009 eV for E || b. Existence of a different 
phonon was indicated by the presence of a third knee at 
E; whose low energy side partner (the fourth knee) was 
not confirmed possibly because of the overlap of free 
carrier absorption. The estimated energy of the second 
phonon, E, — E,, was found to be 0.047 eV for E||aand 
0.022 eV for E|| b. The values of second phonon energy 
are roughly twice the values of first phonon energy. 
This may indicate a two-phonon-assisted pheno- 
menon. The phonon frequencies’ obtained by Raman 
and infrared measurements also support it. 
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Two approximations, namely, the extended random esp ee e Sen ace a par 

(RPA) have been employed to find the direct correlation unction | ne par: 
é ard-Jones and double Yukawa potential functions using the Weeks, Chandler and Andersen 

ace apts Mame ah (USA), 54(1971) 5237] for the division of the pair potential. The ERPA results are found to pore 
exactly identical with the RPA results, indicating that the ERPA does not constitute any significant Bi in a over the ms 
at least for this particular long wave property. Using the compressibility sum rule and through the DCF int € long wave limi 
in the RPA, an equation of state has been formulated by neglecting the small density variation of the WCA diameters during 
integration of the compressibility relation with respect to density at constant temperature. Through this equation of state we 
derive the thermal expansion coefficient, isothermal compressibility, pressure variation of the bulk modulus, specific heat, 
sound velocity of the liquids and also their temperature variations. The agreement between theory and experiment in the 
evaluation of many of these properties shows that this approach forms a very important method in the formulation of equation 


of state. 


1 Introduction 

In a recent work, Evans and Sluckin' have used the 
random phase approximation (RPA) in their study of 
the long wave equilibrium behaviour of simple liquids 
characterized by a Lennard-Jones (LJ) potential. In 
the present work, the author intends to calculate the 
direct correlation function (DCF) in the long wave 
limit for the simple liquids modelled by a L-J potential 
using the extended random phase approximation 
(ERPA) of Evans? together with the Weeks, Chandler 
and Andersen (WCA) prescription? for the division of 
the pair potential. The author extends his calculations 
to include a double Yukawa (DY) potential with 
parameters adjusted to fit a L-J potential. Finally, the 
thermodynamic properties obtained from the RPA 
compressibility equation of state are calculated and 
compared with the experimental results. 


2 Development of Theory 
2.1 Inverse Compressibility in ERPA 
First a WCA division of the Pair potential is made. In 


the WCA prescription the Pair potential g(r) is divided 
into two parts 


P(r) = Polr) + —,(r) ae | 


where Polr) is the reference Part and @,(r) the 
perturbation potential. g(r) includes all the rapidly 


varying repulsive forces and %,(r) all the smoother 
attractive tail. 


Polr) and @,(r) are given by 
Pol) = Ar) — Or.) = G(r) +e 
=@ 


r<Tp 
she wwe) 


280 


9, (r)=9(r,,)= —€ 
= P(r) 


P<T. 
rt ae 


—€éand r,, are the values of the potential minima and 
the separation distance at which ¢(r) is minimum, i.e. 


Oe(r) 
OF spas 


=(): —€&=((r,,) . (4) 

The calculation of the isothermal compressibility 
which is a long wave property of simple liquids is based 
on ERPA for the direct correlation function C(r). By 
treating DCF as a functional derivative of the free 
energy, Evans* showed that in the ‘local density’ 
approximation of the pairwise distribution function, 
DCF is given by 


C(r) = Co(r) — Be, (r)g,(r) . (5) 
Where C(r) and gy(r) are the DCF and the pair 
distribution function respectively of the repulsive 
reference fluid of interaction potential @o(r). 

In order to calculate Co(r), following Jacobs and 
Andersen*, the author approximates it by the sum of 
Cur), the DCF of the hard sphere (HS) fluid at the 


same temperature and density, and the Blip function 
Bir) as 


Colr) = Cudr) + B(r) 


.. (6) 
where 
Bir) = YAr){exp[ —Beo(r)] —exp[—Bedr)}} — ... (7) 
Ydr)=gdryexplBedr)] 


.. a 
where B = 1/(kT), k being the Boltzmann constant and a 


ve i 
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T the absolute temperature. The subscript ‘d’ refers to 
hard sphere quantities where the hard sphere diameter 
is d. It is determined at each temperature and density 
by the WCA requirement, 


[Bindr=o .. (9) 


The compressibility equation of state is then given 
by 

‘d 

I — pC(0)= l ~4np | r?Cydr)dr 


0 


+4npp | r?@,(r)go(r)dr .. (10) 


d 
where p is the number density, C(0) the long wave limit 
of the Fourier transform of C(r). The first two terms on 
the R.H.S. of Eq. (10) is just the inverse isothermal 
compressibility of the HS system. Using Carnahan and 
Starling’ equation of state for HS, it is given by 


(1 + 2n)? +0°(n—4) 
(1—n)* 
where n(=pd?/6) is the packing fraction. 


In the WCA theory’, the reference system radial 
distribution function is given by 


golr)® YAr)expl — Bg (r)] + (1.2) 


Using Eqs (9), (12) and the fact that Y,(r)=g,(r) for r 
>d, the last integral in Eq. (10) reduces to® 


d 
149 | r?Cydr)dr= E88) 
0 


1=4n96 | 9, (r) gdr)r? dr «~ €i3) 
d 

For a Lennard-Jones fluid in the Percus-Yevick (PY) 
approximation’, this integral has been shown to be 
given by® 


12 __yi2 
1=96n.p4 


i- me (4) 


Ms & 


The integrals I{ and /} are aie in Ref. 8 and need 


not be Seeadhsced again. 
It is interesting to note that one can obtain the value 


of I directly from the work of Thiele® 


= -[ (i, tte)? ar 


.. (16) 


Alt ee Ait on approximate result, it gives 
en 


2.2 Compressibility for a Double Yukawa Fluid 

Recently, Foiles and Aschroft'® have shown that a 
double Yukawa (DY) potential provides a good 
physical representation for pair potentials in molecular 
fluids. They have shown that by judicious selection of 
parameters this potential can have very nearly the 
same form as Lennard-Jones potentials in the 
physically important range of separations and can 
therefore lead to essentially the same associated 
thermodynamic properties. The DY potential is given 
by 


E 
Povlr) =~ rexpL — a(x — 1)] —exp[ — b(x — 1)]j eweh 1) 


where x=r/s 


Foiles and Aschroft found that the set of parameters 
E=2.0199¢, s=o, a= 14.735 and b=2.6793 may be a 
good choice for the matching of the DY and LJ 
potentials. Recently'', the author has derived closed 
analytical expressions for the free energy and the virial 
compressibility factor of a DY fluid using the WCA 
perturbation theory*. He now calculates the 
isothermal compressibility of this system in the ERPA 
simply by finding the integral J for the DY case. 

In the PY approximation, J is given by"’ 


See ee . (18) 
where 

1, =24n,, Bel? ;, {19) 
Ke a 24 w= e@G(ac,,) —e°G(be,.)] . (20) 
in nc 24nwBE F(a) —F(b)]) .. (21) 
Bey 


The integrals /,, are given explicitly as: 


1, =4npBe [jet ita ne Par 
= 4npBE {, [setae ont oe 


= 4npBE f kee —a(r/o —1)]—exp[—b (r/o — mt 


r/o 
x oe nN pra 
d d o'iw 


In Eq. (20), 2) is the Laplace transform of 
"Sa Pld, Nw) and is given by Wertheim!? as 


zL(z) 


G(2)= 99, [Li + Se) 


sae (22) 
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Eqs (10) and (27) can thus be used to obtain the 
compressibility values in the ERPA and RPA schemes 
respectively. Values of 1 —pC(0) for the simple liquid 
... (24) argon at different temperatures are given In Table 1. 
Since the LJ and DY potential functions yield more or 


L(z) =12n,[(1 +n,/2)z +(1 + 2n,)] ae toa) 


S(z)=(1 —n,,)?z? + 6n, (1 — Ny)z? 
+ 18;,2z — 12n,,(1 + 2.) 


we 


F(a) and F(b) are given by” less the same values, the author has shown the results 
A ; } from the LJ potential model only. It is interesting to 
DSi Al at iinet note that both the approximations (ERPA and RPA) 
st yield almost exactly identical results for the inverse 
< (1 + sx, expel —€yw)] compressibility. This shows that the use of ERPA does 

ie ae not constitute any significant improvement over the 


I RPA, at least for this particular long wave property. 
(sn +2 expat ~ Xmew) The reason for this is that the integrals { g,(r)r? dr in 
~,(q=0) of Eq. (28) and | g,(r)gdr)r? dr of Eq. (13) 
Bee 2 2 : : ow 
+—j <1+— +——, pexp[2(1 —c,)] differ only by a few per cent. The author has thus led to 
2ewl | 2w (Ze) suggest that the following equation 
2x 2 2 3 
—<xz+—"™+—_, —XmC =e l+n+n*— 
pares + ceppeselal Xmé 3] Came eee eae .. 30) 
pkT (1 —n) 
where 
A=to~ %, +03/2, B=0,—a, and C=a/2  ... (26) obtained by integrating the RPA inverse compre- 
%, %, and &, are given in Ref. 8. ssibility relation (27) can be taken as a convenient 


The use of the DY potential model thus renders the equation of state for the simple liquids, provided that 
evaluation of the integral / in an analytically exactform the diameter d is determined by the WCA criterion 
unlike the LJ potential model where the integrals /" (since at fixed temperature d varies weakly* with 
have been approximated®. density p, hence during integration with respect to p we 
23 Compressibility in the Random Phase Approximation may assume d to be constant). Nevertheless we should 

The RPA results for the compressibility can be remember that the thermodynamic functions calcu- 
obtained from the work of Evans and Sluckin'. The lated through this route do not depend on the ‘softness’ 


result is — of ~0(r), since these involve C o(0) which are hard sphere 
1—pco ="! + 2n) EMD) oigitg ~0) Q7 like by construction (Ref. 9). The RPA critical 
(1—n) constants have already been discussed by Evans and 


where q is the wavevector and Q 1(q =0) is the long wave Sluckin’ and the author need not discuss them again. 
limit of the Fourier transform of ~,(r) 
This is given by! 2.4 Evaluation of Thermodynamic Properties 


01(q=0)= — 15,797¢03 .. (28) The equation of state (30) will now be utilized to 
verify the results for some real fluids, especially with 
; respect to some easily measurable physical properties 
%1(q=0)= — 15.793¢03 ... (29) Of simple liquids such as argon, nitrogen, oxygen and 
methane. The Properties chosen are: (i) isothermal 


for the LJ case and 


for the DY case!!. Actual expression for (q) for the 


OF Gnd sit aiven an compressibility B;, (ii) thermal expansion coefficient 
a 27, (ili) pressure variation of the bulk modulus C, , (iv) 
— “ a 3] 4'mCOS qr, — sin ara | 
|,4%2.0199f expl — alra/o_ » Table !—Comparison of 1—pC(0) Values for Argon 
ae rT mee COS Grp Obtained Using ERPA and RPA 
Temp. 1-eC(0) 
2-2\-1 a es 
+asin ara( +] | ra 
, ERPA RPA 


23.2419 23.2918 
17.9450 17,9075 
17.4988 17.4478 
16.9869 16.9208 
16.6450 16.5676 
16.3865 16.3000 


_ 42.0199 expf — bir,./o — 1)] 
a rT : we To 


+ bsinara( ) |] 


SSBSRE 
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heat capacity C,(/), and (v) sound velocity C in the 
liquid. The expressions for these quantities are as 
follows (in all the calculations, the author has neglected 
the small density variation of the WCA diameter d 
during integrations and differentiations): 


(i) Isothermal compressibility (B,-) 


V 
Br= pel - eC)! uc GY 


where 1 —pC(0) is given by Eq. (27). V is the molar 
volume and R the universal gas constant. 


(ii) Thermal expansion coefficient (a) 
We know that 


mer), "*\ar), 
Differentiating Eq. (30) with respect to temperature 


at constant volume, we get 
a 


HES) 


‘. 
- 
a 


a av | tl —®) 


—C ana ei “| 23) 
+dyoi( 3) . 34) 


“nitrogen and oxygen, and (3) symm 


where 


o 


At w" 1 ¢s a; , er 
d da(1 +513, )+ ody 
” Oo; a, ae 0; N 


Sh ees. er) ML 
aa Ope)? * KBs) 


eee 


... (38) 


a, and d’ can be obtained from Eqs (33) and (34) 
respectively. 

The contribution of the bracketted portion is, 
however, less than | % and can therefore be neglected. 

(v) Sound velocity (C) 

Sound velocity C is given by 
=. oe 

MB; 
where M is the molecular weight of the liquid, and 
7(C,/C,) the ratio of the heat capacities at constant 


pressure and constant volume. C, may be obtained 
from the thermodynamic relationship 


C,—C,=a2TV/By .. (40) 


when B;, %7 and C,(!) are already known from Eqs (31), 
(33) and (36) respectively. 


C2 


... (39) 


3 Results and Discussion 
_ The calculations of the author are reported for (1) _ 


ric 
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Table 3—Variation of Thermodynamic 


P 
Liquid Temp. Pr . 
K 10°*atm™' he 
Cale. Obs Cale. Obs. 
Ar 84 1.78 1.93 3.75 4.45 
86 1.80 2.04 3.69 4.47 
87 1.84 2.10 3.70 4.49 
87.29" 1.84 2.10 3.69 4.50 
88 1.88 2.16 LH Ds 4.51 
89 1.91 P92: 3.71 4.53 
90 1.93 2.28 3.70 4.56 
N, 65 1.88 2.06 4.07 5.21 
70 2.15 2.54 4.12 5.34 
75 2.46 3.03 4.20 5.49 
77.19* 2ST 3.30 4.20 5.56 
80 2.82 3.57 4.30 5.64 
85 3.23 4.11 4.43 5.81 
90 3.68 4.70 4.56 5.99 
O, 65 0.91 1.06 3.19 3.68 
70 0.99 1.19 3:15 3378 
80 1.20 1.49 Sh 3.84 
90 1.46 1.79 3.14 3.95 
90.13* 1.46 1.80 315 3.95 
CH, 102.6 1.88 1.32 3.09 seit 
108.1 2.08 1.46 3.13 3.25 
kEL7* 2.15 1.66 3.10 3.34 
113.7 pI9.9 1.66 £3 Ui 3.40 
119.3 2.42 TT 3.15 3.56 
124.8 2.67 = 3.20 3.86 


* Boiling point of the liquid 


Properties of Simple Fluids with Temperature 


C,(l) Cc C, 
parameter 
cal/mol/K m/sec 
Paes Obs. Cale. Obs Cale. 
9 38 10.49 875.1 860.2 7.92 
9 30 10.51 871.7 846.8 7.84 
9.29 10.52 864.9 840.0 7.81 
9.28 10.52 864.2 — 7.80 
9.28 10.52 856.5 833.1 7.78 
9.26 10.54 851.5 826.4 7.74 
9.24 10.54 848.2 — 7.71 
11.41 13.32 1005.3 975.0 8.20 
11.30 13.44 953.6 923.5 7.98 
11-22 13.50 905.1 873.0 7.79 
11.18 13.56 889.7 —— 7.71 
11.18 13.64 859.4 822.0 7.62 
11.16 13.70 816.4 772.0 7.46 
8 (Is 13.84 778.5 721.0 7.31 
11.63 12.70 1193.9 1104.5 8.66 
11.42 12.77 1146.8 1066.0 8.42 
11.11 12.86 1059.2 988.0 8.01 
10.91 12.99 981.3 909.0 7.66 
10.91 12.98 979.4 —_ 7.66 
12.40 13.1 1386.8 7.95 
12.35 13.2 1333.3 7.81 
12.28 13.45 1315.8 — 7.71 
12.26 13.3 1301.1 7.66 
12.22 13.5 1258.3 7.52 
12:21 13.6 1214.3 7.40 


~~ 


ture'*"'®. Even though the C, parameter has a general 
value of around 8, it may be regarded as a constant 
characteristic of the liquid as confirmed by the 
calculations of the author. 

It is seen from Table 3 that the magnitudes of the 
calculated values of B,, a, and C are more or less in 
good agreement with the observed values and are far 
better than the corresponding hard sphere results of 
Yoshim'*. However, the agreement can be made to 
improve in the cases of a; and C,(/) if the experimental 
values of B; and a, are inserted in Eqs (33) and (36) 
respectively. Another important thing to notice is that 
the trend of variation of the calculated Properties with 


va cases. The 
compressibility B; and sound velocity C reproduce the 


reproduced b 
calculations of the author. Also just like hard Bonn 


fluids'*, the calculated values of Cl) decrease with 
284 


increase of temperature, contrary to the experimental 
results. 

The disagreements may be due to inherent weakness 
of the derivation of the equation of state through the 
RPA. As mentioned earlier, the RPA (as well as the 
ERPA) does not correct the core DCF of hard sphere 
potential due to perturbation?°, which means that 
modifying the form of 9, (r) for r inside the hard core of 
the potential does not change C,,(r). Very recently 
McLaughlin and Young?! have calculated the 
compressibility of argon and neon at several 
temperatures using the RPA and a more sophisticated 
mean density approximation (MDA). However, the 
MDA results are not much different from the RPA 
results. Thus it is gratifying to conclude that the present 
Step is better than the pure hard sphere results and 
forms an important method in the formulation of the 
equation of state for simple fluids. 
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Density, temperature and size of the plasma generated in a plasma focus device, operated = a see parks de 
slow capacitor bank, were determined by measuring the intensity of bremsstrahlung emitted by this p asma - eae “0 
regions. The quantities thus obtained are compared with plasma properties of other plasma focus devices obtained on bigger 
and faster capacitor banks elsewhere, and scaling of plasma characteristics with bank parameters is inferred. 


1 Introduction 

The plasma focus’ device is a high efficiency pinch 
device which produces plasma focus of a few mm 
diameter and a few cm length, with peak density and 
temperature ~10!° electrons/em? and ~10’ K 
respectively and lasts for ~ 100 ns when a high 
powered electrical pulse of a few us duration, generated 
by a capacitor bank, is discharged through it. 

More than 20 laboratories are currently investigat- 
ing plasma focus devices. The present experiments were 
conducted with a comparatively slower and low energy 
capacitor bank with a view to studying the scaling of 
device properties with electrical parameters of the 
capacitor bank. The characteristics of the capacitor 
bank named ‘HEXA’ (Ref. 2) are 16.8 HF capacity, 40 
nH inductance, charging voltage variable up to SOkV. 
The hydrogen filled plasma focus on which the 
experiments were conducted was of Mather type! with 
11.8 cm long central and coaxial outer electrodes 
having diameters of 2.2 and 7.2 cm respectively. 

The plasma formed in the device was characterized 
by experimentally studying the radiation emitted from 
it. The techniques used are described in the following 
sections : 


2 Techniques for Determining Various 

Plasma Characteristics 
2.1 Plasma Radiations 

As is well known, a high temperature high density 
plasma emits copious amounts of bremsstrahlung 
radiation®. Since the plasma formed in the plasma 
focus device is optically thin, it transmits all the 
radiations generated by it and as the temperature of 
this hydrogen plasma is high, the recombination and 
line radiations are insignificant as compared to free- 
free radiation. The intensity as a funciton of photon 
energy of this free-free bremsstrahlung is given by the 
following relation? (assuming Maxwellian distribution 
of electrons): 
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where / is the intensity of the radiation, E (in eV) its 
energy, 7 the temperature (in eV), n the electron 
density, V the total volume of the “hot” plasma, g,, the 
gaunt’s factor for free-free bremsstrahlung, and A a 
constant. As evident from Eq. (1), the radiation emitted 
from the plasma has information about the dimensions 
(V), density (n), and temperature (7) of the plasma. 
Hence by using appropriate techniques (described in 
the next three sections), these quantities can be 
estimated. 


2.2 Temperature 

Since the radiation in the spectral region E> T is 
most sensitive to T, T can be estimated by measuring 
the spectrum of the emission in this region. In the 
present studies, since T is of the order of a few keV, the 
radiation measured is in the soft X-ray range. 

The spectrum of these soft X-rays was estimated by 
the well known foil transmission technique*. The 
attenuation of the X-rays was measured by six 
aluminium filters (thickness, 20, 50, 75, 100, 150 and 
300 jum). Since the ratio of the attenuations through 
any two thicknesses is a function of temperature only, 
the plasma temperature can be deduced therefrom. Six, 
rather than two, attenuating foils were used so as to 
confirm that the distribution is Maxwellian and the 
impurity atoms are negligible (if temperatures obtained 
from all six attenuation ratios are identical). 

As mentioned earlier, since the dense phase of the 
focus lasts only for ~ 100 ns, the intensity of the soft X- 
rays was measured by a detector with fast response 
ume. A plastic scintillation detector (3 mm thick of NE 
102A plastic mounted on a Philips 58 AVP 
photomultiplier tube) with FWHM of 5 ns for a $$Co 
gamma was used for this purpose. The output of this 


detector was recorded on a fast storage oscilloscope © 


(Tektronix 7834), 


2.3 Dimensions 

‘ In the spectral region E < T, the exponential term is 
unity and the intensity of the emitted radiation is not 
very sensitive to temperature [Eq. (1)] but is more 
sensitive to volume and density. For a plasma focus 
device, this condition is satisfied in the optical region 
(E ~ 2eV) and because of the ease and accuracy with 
which visible light could be measured, this region was 
chosen for measurements. The intensity of light in the 
spectral range E to (E + AE) is given by 


AE 
T= An’V orp 855 pee 


As can be seen from Eq. (2), to estimate n or V from 
the light output, the other parameters should be 
known. Generally, volume is measured by photo- 
graphing the plasma. But since the hot phase of the 
S plasma in the plasma focus device lasts only for ~ 100 
4 ns, determination of the plasma volume would call for 
= an ultrafast camera®. Such cameras are not readily 
available as they are very expensive. A new and simple 


hout the use of a fast camera, for a plasma having an 
; of symmetry, is described in what follows. 

‘he plasma in the plasma focus device is 
imately a cylinder. If this plasma is viewed by a 
nato or (Fig. 1) from direction perpendicular to 
dri cal axis, the intensity viewed by the detector 
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method for deducing plasma dimensions and density 
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same plasma at the same point is viewed through a pin 
hole collimator (with diameters much less than plasma 
radius), again perpendicular to the axis of the 
cylindrical plasma, aligned in such a way so as to view 
only the central part of the cylinder, then the intensity 
(1) will be 
2 
In = Ar Bry nr .. (4) 


where a is the area of pin hole collimator. The radius (r) 
of the plasma can then be deduced from the ratio of the 
two intensities given by the following equation: 


r= —— oy 


If the plasma is also viewed by a pin hole collimator 
aligned axially (Fig. 1), the intensity (Jy) is given by 


An* AE 
ly = 855 Fare! ...(6) 


Then effective length (/) of the plasma can be estimated 
from the relation: 


4na lyl 
l= ae : ee f 


Thus the plasma dimensions can be determined. | 

In the experiments conducted by us, the light emitted 
from plasma was again viewed by a fast detector—a __ 
chrome filter (bandwidth, 4500-5500 A) mounted ie ae 
photomultiplier tube (Philips 58 AVP) and the signa & 
were again recorded on the Tektronix fast oscillo: ; 
The slit collimator consisted of a pair of slits of 0 
width, while the pin hole collimator was 
two ies ‘of 0.1 mm diameter. The ali 
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Table 1—Properties of the Trombay Plasma Focus De 
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obtained ones. This suggests that the present th 
about plasma focus phenomenon are not | 
adequate and require further refinement. 
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An apparatus to maintain the equilibrium state of the freezing/melting of tin and zinc over long periods of time and hence 
useful for calibration of platinum resistance thermometers is described. Uniform temperature is established inside the well of 
the graphite crucible and the slope of the flat portion of the melting/freezing curve is +0.2 mK and lasts for nearly 90 min. 


1 Introduction 

The melting/freezing points of pure. metals are 
assigned definite temperature values on_ the 
International! Practical Temperature Scale (IPTS). The 
apparatus described here realizes the thermal 
equilibrium state of the freezing/melting of the spec- 
pure tin and zinc under | normal pressure to a high 
precision (+0.0002°C) required for the calibration of 
standard platinum resistance thermometers (SPRT). 
This work is a part of the programme to realize and 
maintain a National Temperature Scale according to 
the definition of the IPTS.' ~ 


2 Experimental Details 
2.1 Apparatus 

The basic design of the apparatus for realizing mp/fp 
of pure metals is similar to that followed in other 
laboratories* ~°. Fig. | shows the tin or zinc point 
furnace containing the metal holder. 

The graphite crucible, made of pure synthetic 
graphite, was designed to hold spec-pure metal, about 
1500-1700 g, in the annular space and is contained ina 
pyrex glass mould. The graphite crucible also has a 
central re-entrant well for admitting standard 
platinum resistance thermometer through a thin 
walled SS tubing closed at the lower end. Four stainless 
steel radiation shields, acting also as heat shunts, are 
supported by the SS tubing in the region above the 
crucible. The inner space in the pyrex mould below and 
above the graphite crucible is filled with fine glass 
wool. is 

A suitably designed brass cap, containing a water 
condenser, is sealed at the top of the pyrex glass mould 
with piecin wax. The provision of the condenser 
maintains the piecin joint at room temperature 
preventing any leakage. The glass mould is connected 
to pumps and a mercury manometer through a large 
thermostated volume. The hermetically leak-tight 


ate, 


system can be filled with pure dry argon/helium gas. 


The graphite crucible is in the central uniform 
temperature region of the furnace. The metal-holder is 
surrounded by a heavy brass cylinder of large thermal 
inertia. The graphite crucible is suitably insulated 
against thermal losses as shown in Fig. 1. 
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Fig. 1|—Furnace assembly 


289 


INDIAN J PURE & APPL PHYS, VOL 21, MAY 1983 


The power supply (230 V, 50 Hz) from mains is fed 
into an automatic servo voltage regulator, its output of 
230 V+0.5% is next fed into another servo voltage 
regulator and the stabilized output (230 V + 0.257) 
finally to the furnace through a variac. 


2.2 Filling of the Graphite Crucible 

The cleaning of the metal-holder parts and filling of 
the graphite crucible under vacuum with 1500-1700 g 
of pure metal follows the procedure adopted by the 
previous workers* °. 

The system after thorough evacuation and 
degassing is filled with pure dry argon/helium gas at | 
normal atmospheric pressure when the crucible is close 
to the mp/fp of the pure metal. 


2.3 Melting/Freezing Curves of Tin and Zinc 

Starting with an initial temperature slightly lower 
than the melting point, the input power is increased to 
obtain a temperature rise of 0.5°C/hr. The level of the 
molten metal is about 170-180 mm above the top of the 
sensor of the SPRT in the re-entrant well. The design of 
the metal-holder and furnace assembly provide 
uniform temperature over the entire length of the space 
in the re-entrant well within the crucible. 

Measuring the resistance of the SPRT at intervals of 
few minutes, the plot of the resistance versus time is 
drawn and the flat portion of the curve indicates the 
melting process. 

In Fig. 2; V,, V2, V3 and V, indicate decreasing 
voltage inputs, i.e. decreasing raics of heating. The 
corresponding melting curves are illustrated qualit- 
atively by broken curves (m,---m,, M---M,, M3-—-m, 
and m,4-—m,) yielding increasing periods of melting. 
Obviously, when the rate of heating (V,) is rather 
high, superheating and negligible width of flat region 
(m,-—m,) are obtained. 

For realizing the freezing point one starts the 
experiment a few degrees higher than the mp of the 
metal. The molten metal in the crucible is allowed to 
cool 1-2°C/hr. 

_In Fig. 3, the power i put to the furnace is set at 
ee | ‘A eo iA be ip V3, Va, Vs (where V, 
tag AL: corresponding rates of 
cooling of the metal would be faster for V, and slower 
Successive values of V,,. The corresponding freezing 
curves are illustrated qualitatively by broken curves 
A---A, B---B, C---C, D—D, or D’---D’ and E---E, or 


curve B---B lasts for more than an hour corresponding 
to slower rate of cooling. 

For still slower rate of cooling for the voltage setting 
V,, the flat region of the curve C---C may last even less 
than an hour. It would be observed that here the 
freezing was introduced in a greater undercooled state 
than in case of curve B---B. If the freezing was 
introduced at still lower undercooled state, the curve 
traced would be C’---C’, the latent heat then released 
was not sufficient to raise the temperature of the 
assembly to the freezing point of the metal. The 
remaining freezing point curves are similarly explained 
for voltages V, and Vs. 

As tin tends to super-cool by more than 10°C, the 
molten metal is allowed to undercool only by 0.5 to 
1°C before the nucleation is introduced in the melt. 
This is achieved by taking out the SPRT and inserting a 
brass rod (at room temperature) for intimate contact 
for 2 to 3 min with the inner surface of the well and then 
putting the SPRT back. 

Uniform temperature is soon established in the well 
within 10-20 min and usually the flat portion of the 
freezing curve for tin, as monitored by the SPRT, lasts 
for 45 to 90 min (Fig. 4). 

In the case of zinc, the super-cooling effect is less 
than a degree celcius, so introducing the brass rod for a 
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Fig. 2—Melting point curves 
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Fig. 5—Typical freezing point curve for spec-pure zinc 


few seconds was adequate to induce nucleation. Fig. 5 

shows a typical curve at the zinc point for one of our 
SPRTs. 

- The inset curves A’ and B’ in Figs 4 and 5 show the 

flat region of the curves corresponding to the freezing 

process on an expanded scale. Any fluctuation of the 

’ resistance of the thermometer in the flat region is small 

and less than 0.0002°C in the case of both tin and zinc. 


2.4 Resistance Measurement 
The resistance in these experiments was measured 
with an ER thermometer bridge (Leeds & Northrup 
~ Co, USA) and a Null-detector (Leeds & Northrup Co, 
USA 9828-2) to a precision of 10~° Q. Usual 
corrections for bridge errors and self-heating effect 


were applied. 


than a few mm of mercury by controlling 
ithin +2°C. In our system, the 
to same ambient room 

nained constant for several months. 


With helium, the pressure drop noticed was less than 2 
mm in 4 weeks. 

During the flat plateau of the freezing curve, the 
ambient temperature around the aluminium head of 
the SPRT was raised by about 2°C. It was observed 
that this had negligible effect on the resistance of the 
thermometer. This effect in any case did not exceed 0.2 
mK showing the height of the melt (about 220 mm) in 
the crucible was sufficient to nullify the cold- 
conduction down the stem, leads and radiation effect 
from the ambient temperatures. 

In some of the experiments, we induced 
freezing/nucleation in the tin melt by using a short 
brass rod, about 50-60 mm instead of brass rod of 220 
length, so that freezing process was initiated in the 
melt, close to the well only in the bottom portion or 
near the top region. This was done for the introduction 
of non-uniform freezing in the melt on its surface in the 
proximity of the re-entrant well. Contrary to 
expectation, no larger slope or fluctuation in 
temperature in the flat region of the freezing curve was 
observed. Perhaps the re-introduction of SPRT itself 
(with stem usually at lower temperature) helped rapid 
establishment of the uniform temperature in the well 
and freezing progressed uniformly. 

Only small variations in temperature of the order of 
0.2 mK in about 45-90 min during freezing process is 
observed. This also indicates that the metal (source— 
Johnson Mathey & Co, UK or Commin Co, Canada) 
was highly spec-pure. 

It has been observed that for realizing a good flat 
plateau during melting, the rate of heating has to be 
much slower than the rate of cooling required to realize 
a good plateau during freezing. Apparently the reason 
may be that the melting process of the metal begins on 
the outer surface and may not be uniform over the 
entire height. However, this could not be investigated. 


4 Conclusion 

The main feature of the apparatus is the 
incorporation of a water-condenser in the cap of the 
metal-holder. This helps hold the piecin wax joint leak 
tight and the ‘Fixed Point’ functional continuously, 
over a long period. 

The effects of the varying rates of heating, cooling 
and of undercooling before nucleation on the flat 
plateau of the mp or fp curves discussed here show that 
it is necessary to have optimum control on the rate of 
heating or cooling for getting a good plateau. 

The experimental set-up is good enough to allow 
realizing the melting/freezing point for a duration of 
45-90 min and the slope of the plateau of the curve not 
exceeding 0.2 mK. 

Due to varying oxidation state of the PRT sensor 
affecting the stability of a SPRT, the reproducibility of 


291 


INDIAN J PURE & APPL PHYS, VOL 21, MAY 1983 


the melting/freezing points of pure metals can be relied 
to 0.5 mK only. 
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Five matrices of different ionization potentials, containing traces of different elements were prepared. The optimal arcing 
conditions were determined. Application of the magnetic field revealed a relation between the ionization potential of the matrix 
and the magnetic enhancement factor of spectral lines. Matrices of high ionization potentials caused remarkable magnetic 
enhancement of the spectral lines, while those having low ionization potentials showed decrease in the intensity of spectral lines 
On application of the magnetic field to the arc. Moreover, it was found that the changes in the temperature and the electron 
Pressure of the arc plasma, induced by the application of the magnetic field, may be one of the reasons for the magnetic 


enhancement of spectral lines. 


1 Introduction 

The application of external magnetic fields of 
different configurations to the dc arc, for increasing the 
spectrochemical sensitivity, has attracted the attention 
of many investigators’ *. In many cases, the magnetic 
field enhances the emission of spectral lines and hence 
provides new possibilities for the detection of 
extremely small amounts of elements. Vukanovic et 
al.'- attributed the observed intensity enhancement to 
the increase of the residence time of the particles in the 
plasma. 

Leushacke and Nickel® used trace amounts of 
elements in a graphite base and a homogeneous 
magnetic field to study the magnetic enhancement of 
spectral lines. They found that the maximum 
intensification was obtained when the concentrations 
of the elements were near the limit of detection. 
Harisonov and Zadgorska® studied the effect of a non- 
homogeneous magnetic field on an arc burning in 
different gases differing in their ionization potentials. 
They found that the enhancement of the spectral lines 
by the application of the magnetic field depends on the 
ionization potential of the gas and attributed it to the 
increase of the number of ions in the arc plasma. 

As an extension in this field, it is of interest to study 
the effect of matrices of different ionization potentials 
on the magnetic intensification of spectral lines. 


2 Experimental Details and Results 
The arc was placed in a radially symmetric non- 
homogeneous magnetic field. The field was produced 
bya water-cooled electromagnet coil provided with an 
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iron core whose upper end has a conical shape of solid 
angle z/2. The axis of the arc was in coincidence with 
that of the magnet. The magnetic induction B and its 
vertical gradient dB/dz, which could be varied by 
changing the current in the magnet coil, were measured 
at different heights from the pole. 

A medium quartz spectrograph Q24 (Carl-Zeiss 
Jena) was used, with a slit width of 10 xm, to obtain the 
spectra. The spectra were recorded on ORWO 
photographic plates type WO 3 and the intensities of 
the spectral lines were evaluated from characteristic 
curves obtained by the preliminary curve method. 


2.2 The Elements under Investigation 

The elements selected for the present investigation 
can be divided according to their volatility’ into three 
groups: 

(a) volatile: In, Sn, Ag, Cu and Zn, 

(b) medium volatile: Al, Mn, and Mg, and 

(c) involatile: B and Be. 

In each group, representative elements with different 
ionization potentials are chosen. A mixture of these 
elements is added to spec-pure graphite powder in such 
a concentration that for each element the most 
sensitive lines appear in the spectrum of the arc. This 
mixture will be referred to as the ‘O’ mixture. Table 1 
gives the concentrations of the elements in the ‘O’ 
mixture, their ionization potentials and the 
wavelengths of the spectral lines used for the 
investigations. 


2.3 Factots Affecting the Magnetic Enhancement 
of Spectral Lines 
The magnetic enhancement of spectral lines is a 
function of some experimental variables; these are: the - 
concentration of the element, the arc current, the shape 
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of the electrodes, the arc gap, and the applied magnetic 
induction. In order to achieve the optimal working 
conditions, the effect of all these variables on the 
enhancement factor I/I, was studied. The variables / 
and J, are the line intensities corrected for background, 
with and without magnetic field respectively. This 
investigation was carried out using a graphite matrix. 
The Snl line at 2839.99 A was taken as a test line. The 
results of this investigation can be summarized as 
follows: 

2.3.1. Effect of the concentration of the element— 
The previously mentioned ‘O’ mixture of elements was 
diluted with spec-pure graphite powder in different 
ratios and then //I, of the test line was determined. The 
result is shown in Fig. 1, which shows that the 
maximum enhancement corresponds to a dilution of 
1:9. This result is in agreement with that reported by 
Leushacke and Nickel’. 

2.3.2. Effect of arc current—The effect of the arc 
current on the enhancement factor I/I of the test line is 
shown in Fig. 2. It is obvious that the optimal value of 
the arc current is 8 A. 


Table 1—TIonization Potentials Spectral Lines and 
Concentrations of the Elements in the ‘O’ Mixture 
Element lonization Spectral Conc. 
potential line, A fi 
eV 
In 5.79 3039.36 0.0012 
Sn 7.34 2839.99 0.0053 
Ag 7.57 3382.89 0.0007 
Cu 7.72 3273.46 0.0007 
Zn 9.39 3345.02 0.0012 
Al 5.98 3092.71 0.0032 
Mn 7.43 2798.27 0.0050 
Mg 7.64 2852.13 0.0006 
B 8.30 2497.73 0.0036 
Be 9.32 3321.34 0.0022 


ss 
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0.5 
0 y] 4 6 8 10 12 14 
Dilution ratio 
Fig. |—Dependence of the magnetic enhancement 
factor I/Io of the test line SnI 2839.99 A on the dilution 
ratio of the ‘O" mixture [Arc current: 6A, crater 
diameter: 4mm, crater depth: 4mm, cathode: 
pointed, are gap: 2 mm, magnetic induction : 11 mT, 
height of arc gap from the magnet pole: 5 cm] 
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Fig. 2—Effect of arc current on the magnetic 
enhancement factor I/I9 of the test line SnI 2839.99 A 


Table 2— Dependence of the Magnetic Enhancement Factor 
I/Iy of the Test Line Sn] 2839.99A On‘the Shape of the 
Anode 


[arc current: 8A, cathode: pointed, arc gap: 2 mm, magnetic 
induction: 11 mT, height of arc gap from the magnet pole: 5 cm] 


Crater I/I, for crater diameters of 
depth, —— 
mm 2mm 4mm 5.2mm 
zZ 1.30 1.60 1.71 
4 1.88 2.55 2.65 
6 1.50 2.11 2.21 


EEE 


2.3.3. Effect of the shape of the electrodes—The 
effects of the lactors governing the shape of the lower 
electrode (anode) namely, the internal diameter and the 
depth of its cavity on the enhancement factor I/I, of the 
test line were studied. The results are given in Table 2. 
From these results, one can see that I/Ig exhibits a 
maximum at a crater depth of 4.0 mm which increases 
gradually with the increase of the internal diameter of 
the crater. Meanwhile, with increase in the crater 
diameter it was found that the arc started to wander 
and consequently fluctuations in the spectral line 
intensity were observed. For the purpose of the present 
investigation, the optimal crater diameter was taken to 
be 4.0 mm. 

Three types of counter electrodes were used. The first 
two are the preformed 100 u and L3963 electrodes. The 
third electrode is a graphite rod 5 mm diameter pointed 
at its end to form a cone with solid angle of 24°. The 
effect of the shape of these types of counter electrodes 
on I/Io of the test line and the corresponding relative 
Standard deviation, using five sates 
pin ar tive Measurements at an 

Pp arc current of 8A indicated that the highest 
enhancement with the least Standard deviation j 
obtained using the pointed counter electrode. a 
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Fig. 3 shows the effect of the arc gap on the 
enhancement tactor //I, of the test line. It can be seen 


that maximum enhancement is Obtained with an arc 
gap of 4.0mm. 


2.3.4. Effect of the magnetic field on the intensity of 


Spectral lines—Fig. 4 shows the enhancement of the 
intensity of the test line as a function of the magnetic 
induction B at different heights z from the magnetic 
pole. It can be seen that the maximum intensification of 
the spectral lines is achieved at a magnetic induction of 
7 mT and a height of arc gap from the magnet pole of 5 
cm. This corresponds to an induction vertical gradient 
of 1.9mT/cm. It was found that this result holds true 
for the spectral lines of all elements under investigation. 


24 Effect of Matrix on the Magnetic Enhancement of Spectral 
Lines and on the Changes of the Plasma Parameters 

The matrices chosen for this investigation are BaO, 
CuO, SiO,, As,O, and graphite. The choice is based 
on the ionization potential of their cations, which are 
5.21, 7.70, 8.35, 10.05 eV respectively and that for 
graphite which is.11.30 eV. The charge of the electrode 
consisted of 70% matrix, 20% spec-pure graphite 
powder and 10% of the ‘O’ mixture. For each matrix, 
eight electrodes were prepared; four of them to be arced 
with magnetic field and four without. So, each intensity 
value was obtained as the average of four readings. The 
optimal arcing conditions, achieved previously, for the 
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Fig. 3—Effect of arc gap on the magnetic 


enhancement factor I/Io of the test line SnI 2339.99 A 


12 % 
induction (mT ) 


Magnetic 


Fig. 4—Dependence of the magnetic enhancement 

factor I/I, of the test line Snl 2839.99A on the 

magnetic induction B and the height z of the arc gap 
from the magnet pole 


graphite base were conducted to study the effect of the 
different matrices on the magnetic enhancement of the 
spectral lines of the different elements as well as on the 
changes in the plasma parameters of the arc, induced 
by the magnetic field. The magnetic enhancement 
factor I/Io, for each element in the mixture, was 
determined from the spectra. Moreover, the excitation 
temperature and the electron pressure with and 
without magnetic field were measured. The excitation 
temperature was determined from the intensity ratio of 
the zinc atom line pair Znl 3072.06 A/ZnI 3075.90A 
[Ref. 8]. The electron pressure was determined from 
the intensity ratio of the ion-atom line pair of 
magnesium MgllI 2795.53 A/Mgl 2779. 83A [Ref. 9]. 
The results of this investigation are shown in Tables 
3 and 4. Table 3 gives the effect of the matrix on the //Io 
ratio for the different elements. The ionization 
potential of the matrix is also given (Table 3). Table 4 


Table 3—Effect of Matrix On the Enhancement Factor //I) of Spectral Lines of Different Elements 


I/Io 
Volatile Medium volatile involatile 
OI fest ood shield 
Matrix 7 a Sn As Cu Zn Al Mn Mg B Be 
0.36 0.44 0.95 0.34 0.63 0.49 0.58 0.26 0.54 
“ae nea re 0.40 0.40 a 0.37 0.73 0.60 0.67 0.49 n.d. 
= yh 0.54 0.22 0.19 0.61 0.81 0.70 0.70 0.51 n.d. 
Si Bt has 146 1.01 1.01 0.70 1.49 1.81 1.67 0.73 0.78 
eae? pe 431 3.80 2.31 2.11 1.00 4.22 2.00 3.11 1.31 1.30 
: ‘d. means the line was not detected in this matrix. 
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Table 4—Effect of Magnetic Field on Temperature and Electron Pressure of the 
Arc Plasma using Different Matrices 


T,K P1072 
Matrix I/Io a 
$n12839.99A 0 7mT AT 0 7mT AP, 

BaO 0.36 5620 5760 140 2.03 1.22 —0.81: 
CuO 0.40 6100 6260 160 3.06 2.34 —0.72 
SiO, 0.54 6120 6220 100 0.99 0.40 —0.59 
As,0, 1.46 6100 6200 100 0.92 1.21 +0.31 
graphite 3.80 6210 6510 300 1.08 2.48 + 1.30 


6.2? ee ee eee 
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Fig. 6—Dependence of the magnetic 
factor I/Io of the test Ii 
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dependence on the ionization potential of the matrix. 
In the presence of matrices of low ionization potential, 
the application of the magnetic field causes a decrease 
in the electron pressure. These matrices depress the 
intensity of the spectral lines of all elements present. On 
the other hand, matrices of high ionization potential 
cause an increase in the electron pressure when the 
magnetic field is applied. These matrices enhance the 
emission of the spectral lines. A relation is obtained 
between the relative change in intensity Al/I of the test 
line SnI 2839.99A and the change in the electron 
pressure AP,, induced by the magnetic field. This is 
illustrated in Fig. 7. 

From the results of this investigation, it may be 
concluded that one of the reasons for the magnetic 
enhancement of spectral lines can be the changes in the 
plasma parameters induced by the magnetic field. 
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A microprocessor-based system for internal friction measurements is described. The system is built around a National 
Semiconductors’ SC/MP microprocessor and can handle set-ups covering a wide range of frequency and temperature. 
Experimental parameters like initial temperature, temperature increment/decrement, number of temperature steps, etc. can be 
specified by the user and thereafter the experiment is fully controlled by the system. 


1 Introduction experiments covering a large number of relaxations 
The study of internal friction in solids provides a covering point defects, dislocations and grain 
convenient means of obtaining information about the boundaries. Similar automated systems have been 
thermally-activated defect motion in solids. Examples _ described earlier? ~*, but the present one differs from 
are the Snoek, Zener and Bordoni relaxations which these in that the logic is not hardwired. Functional 
occur due to interstitials, substitutionals and requirements are met through software programming, 
dislocations respectively’. The type and orientation of _ thus lending flexibility and ease in system design. 
the defect can be identified by the characteristic ; . 
dependence of the internal friction Q~! (or damping) 2 Block Diagram and Control Algorithm 
on temperature and frequency, given by 
Q°* = Aat/(1 + w?1?) rod) 
where is the angular frequency of oscillation of the 
sample and t the relaxation time for the process giving 
rise to internal friction. A is a constant which contains 
information regarding the concentration of defects, 
their strength to cause damping, etc. 
The relaxation time Tt is Strongly temperature 


dependent through the Arrhenius relation voltage is directly fed to the CPU from the multiplexer — 
eee!) | ---(@) using ADC-1. The CPU, in addition, has main memory _ 
where H is the activation energy for the process, T the (where the program is stored) and a keyboard 
temperature, kg the Boltzmann constant, and to the printer as the 1/O devices. 
frequency factor. It may be seen from Eq. (1) that,for Before proceeding with the description of the control 
€ach process, the maximum (zt = 1) will occur for a algorithm, we present the sequence of operati aes 
Particular choice of frequency and temperature range performed at each temperature. These ar nee 
[as dictated by Eq. (2)]. The wide variation in w and T | . Peri 
entails the use of different experimental arrangements. 
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Controlling the temperature, (b) exciting oscillations to 
3 Shar tn a preset level, (c) recording t cok ane 
Further, to obtain the Debye curve completely [Eq. of the Prin a the a PTE ee ce 
(1)], one has to perform the measurement of damping period, and finall Ms tions, (d) x ring 
Saletan cae be, com extreanely tedious Printing out : values of a 
especially if one is dealing with low frequencies. . ne period. The sequence (b) tc 
large amount of data that has to be handled at Pr Or cet uae thrice at each temperature b 
underscores the need for some kind of automation. The pare 
This paper describes a mic based ° 
measurement system that facilitates : 
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2 _ OD 
TRISTATE BUFFER TRI-STATE BUFFER 
onvtollP 12 BIT ADC 2 


12BIT ADC 1 Lend gf 


start conv 


PEAK DETECTOR 


from cap transducer 


from thermocouples 


TO DECODERS & MEMORY 


RAM EPROM 
512x8 256 x8 


O/P PORT 1 O/P PORT 2 


RI-STATE BUFFER RI-STATE BUFFER 
to temp. controller to excite 


end of conv. 


to interrupt of UP 


Fig. 1—Microprocessor-based internal friction measurement system: Block diagram 


other tasks like excitation, data acquiring, com- 
putation, etc. are subroutines incorporated into the 
interrupt service routine and are executed only when 

need 

The flow chart of the temperature control routine is 

shown in Fig. 2. Upon ‘cold starting’, the system 
initializes its registers and waits for the user to enter the 
experimental parameters. These are the starting 
temperature, temperature step size, number of 
temperature steps, sensor type, etc. (however as a 
default option, certain typical values for these 
parameters are assumed). Once this step is completed, 

the system latches each of the four temperature 
channels (inputs) sequentially, and at each channel it 
initiates analog to digital conversion of the amplified 
sensor voltage. These digital values are then stored. All 

these voltage readings are converted into temperature 
(linearized) and a difference table of A7’s (difference 
between current temperature at which a reading is to 

be taken, Te: and Tactuai) is prepared. The sign of AT is 

used to switch ON or OFF the power to the heater in 

| each channel. This process is continued till any one of 
the AT7’s is zero. This implies that the particular set-up 
has reached the set temperature and is ready. Next the 
DATA flag is set and the interrupt of the 
r is ‘enabled’ (The data flag has been 
sed so that subsequently when AT = 0 for the other 
channels, the microprocessor does not attend to them 
Hi is disposed of). The appropriate 
i given an excite input and 
tput) is 


& 


ENTER EXPTL. PARAMETERS 
INITIAL TEMP, TEMP STEP NO. OF STEPS 


READ TEMP SENSOR VOLTAGE 


MP-VOLTAGE 
ACTUAL 


LINEARIZE T 
TO GET 


TURN ON HEATER 


TURN OFF HEATER 


(S DATA FLAG SET? 
SET DATA FLAG AND 
ENABLE INTERRUPT 


control routine and main program of the 
system: Flow chart 


Fig. 2—Temperature 
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latched. The system then continues with the 
temperature control algorithm. 

It may be noted that the actual collection of data 
takes a short time (+30 sec for a frequency of | Hz) 
compared to the time taken to raise or lower the 
temperature of the sample. This fact is utilized in 
attending to each set-up (‘servicing’). The servicing 1s 1n 
terms of which set-up has reached the desired 
temperature first. However, in the event of two 
channels showing AT = 0, the set-up with the lower 
address designation is given priority. Thus in the long 
run, all experimental set-ups will get equal attention 
provided the thermal time constants are not grossly 
mismatched. 

We now proceed to discuss the sequence of 
operations once a particular channel has reached its 
current 7... A method to ensure smooth excitation of 
oscillations is that the sample be excited repeatedly till 
the p-p amplitude reaches a preset value (This can be 
varied if needed). This is accomplished by the use of 
another flag, the data acquisition flag (DAQ) in the 
interrupt service routine (Fig. 3). A check is first made 
on the DAQ flag and if it is not set, the system proceeds 
with the excitation, i.e. an excite pulse is generated after 
each positive peak. Otherwise amplitudes of 68 
positive and negative peaks (34 p-p values, A,) are then 
recorded, each recording being an interrupt by the 
data. 

After acquiring all the data (ie. the p-p values 
corresponding to the channel which has reached T;.), 
the first four peak values (2 cycles) are ignored as these 
May correspond to transients. The value of the log 
decrement 0 is then calculated from the remaining data 
points by an averaging formula 


= (1/N) in = ‘A= lo ai ...B) 


Next, the time period is read in and the results are 
printed out along with the channel number. The DAQ 
flag is then reset and the experiment is repeated two 
more times. At the end of three iterations, the DATA 
flag is reset. The system finally checks for any 
erroneous condition (e.g. poor vacuum) and prints out 
a condition code before Proceeding to increment 
temperature and perform data acquisition the next 
time it meets AT =0 (either for this channel or any 
other channel). The condition code has been used 
purely as an indicator and no corrective action is taken 
Before setting the next temperature ona particular 
channel, a check is made whether the experiment is 
over or not, i.e. whether the temperature cycle has been 
gone through fully. This is done as follows: The 
number of temperature steps is compared with the 
user-specified value. If the two are unequal, the step 
count is incremented and the experiment proceeds. If 


ED leis 


i=W41 


SAVE REGISTERS 


EXCITATION OVER? 
NO CHECK DAQ FLAG YES 
GIVE EXCITE PULSE ACQUIRE DATA 
DATA 
AQUISITION OVER? = 
YES 


COMPUTE @’ 
so —<Gone 2 tes 
YES 
er) ee aI Grn 
0 ; 
Fig. 3—Interrupt service routine (All RETURNS are to the main 
program) 


they are equal, the temperature increment 57 is 


' “ 
is 


- 


checked to see if it is negative (This option has been 
included so that a heating-cooling type run can be 
done). If it is positive, the step number is set to zero and 
the increment is negated. Otherwise the experiment is 2 
Over and the message is printed out. ; = 
The computation and printout take only about five 3 
seconds so that temperature control is quite effective. _ 
However, higher speeds of data acquisition may not be __ 
compatible with such an interrupt service structure.In 
such cases, a direct memory access (DMA) controller _ 
will have to be added and the microp or will be 
interrupted only after all the data have been acquired. __ 
This entails only minor changes in the interrupt service 
routine. | war ee 
It may be observed that the initialization of 
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drawback has been taken care of to a certain extent by 
having a ‘LOCAL’ mode on each of the temperature 
controllers which will keep the temperature at their 
current values. The fourth set-up can now be connected 
and the system appropriately re-initialized before 
changing the controllers to the ‘LINE’ mode. 


3 Hardware 


3.1 Transducer 


The typical strain amplitude of the sample is in the 
range 10° °-10~*. A capacitive transducer has been 
employed in our set-up and the corresponding change 
in capacitance works out to be in the range 102-107 * 
pF. We have used the ratio transformer principle to 
measure these small changes in capacitance’. The ratio 
transformer and associated circuitry are depicted in 
Fig. 4. The rms current at the point A is given by 


in = (Vo/./20(C, — C)@o .. (4) 
where +V osin@ot is the voltage appearing across 
the secondary of 7,. The output is thus proportional to 
the difference in the capacitance of the two arms (which 
are attached appropriately to the mechanical part). 
Slow variations in (C, — C,) can be tracked provided 
the frequency of such variations is much smaller than 
@o. We have chosen a frequency of 10kHz which 
allows tracking variations up to 100 Hz. 
__- The current i, is typically a few picoamperes for C, 
a Se ~ pF and V)=2V. This current can be 
| 1 into a voltage which is measured using a 
Pr. ampli 1er. As shown in Fig. 4, the current to 
ye conversion is achieved through the FET input 
np. Sean has an input resistance of 10!?Q 
2¢ >A input bias current. This allows for a 1 MQ 
ck resistor to be placed without saturating the 
Ad circuit layout was made quit ae 
atch. 
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microprocessor at 110 baud to simulate the operation 
of a teletypewriter. 

The 16 bit address is formed from the address pins 
(12 LSB address bits) and from the data bus. The 
address decoding is made simple by reserving address 
lines for specific devices, as for example, the address 
line 15 is reserved for the temperature sensor channels 
and 13 for the capacitor transducer channels. 

The input channels have been divided into four 
groups; the thermocouple input channels which are 
slow and require an integrating type ADC and the fast 
channels which are connected to a fast ADC via a peak 
detect/peak locate circuit using ADC’s with conversion 
times of ~ 10 msec. While this at present caters for the 
low frequency experiments, a 12 bit successive 
approximation ADC (conversion time ~ 20 psec) is to 
be used with the medium frequency set-ups. 

The peak locate signal from the peak detector’ is 
used to initiate conversion on the signal channels and 
the same pulse is fed to a timer circuit to measure the 
time period of the oscillations. The conversion on the 
temperature channels is initiated under software 
control. The end of conversion from the signal channel 
ADC is used to interrupt the processor to read in the 
peak data. 

The temperature controllers are of the on-off type 
and the heater powers are switched on or off depending 
on the control signal from the microprocessor. This is 
achieved by means of a triac circuit and the control 
signals for the temperature controller are latched on to 
output port 1. The output port 2 forms the excite signal 
port and, whenever a particular set-up is to be set into 
oscillation, the bit corresponding to that set-up is made. 
high. In addition, a digital input port 3 is provided to ; 
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Fig. 4—Ratio transformer bridge and associated circuitry 


are stored in a look-up table and temperature 1s 
computed by the formula 

T, =((T, — T)V. — V/V 5 — VIN + Ti ey, 
where V ,, V;, T,;, T; are the voltages and temperatures 
at the end points. and 7, is the temperature 
corresponding to the voltage V,. The above scheme 
gives an accuracy of about 0.5% for the K type 
thermocouple. Temperature control is achieved by 
comparing this temperature 7, with the set 
temperature value stored in the RAM. The 
temperature control output is set low or high 
depending on the result of this comparison. 

The value of 0 is calculated by the relation given in 
Eq. (3). The sums are computed using multibyte 
arithmetic and then the logarithm of these sums is 
calculated. The sums are first represented as a fraction, 
X (the sums are first left justified by shifting and the 
number of shifts is noted). The logarithm is computed 
by multiplying the given number by certain constants 
whose logarithms are known, such that the product 
approaches unity, ie. 


Bivtsls..fr=1 ... (6) 
and hence 
InX = —(alnf, +binf, + ...+minf,) wiht) 


The factors f;, f,, etc. are chosen such that the 
multiplication becomes a trivial shift and add. The 
values of their logarithms are stored as decimal 
numbers and the result is obtained directly as a BCD 
number which is printed out after conversion to the 
ASCII code. 

The conversion of numbers from binary to BCD and 
then to ASCII is done very easily by the use of the 
decimal ADD instruction of the microprocessor and 
adding 30 to the resulting BCD code. 


5 Results and Performance 


The reproducibility of the 6 values was found to be 


better than 2%. This was determined b 
0 y repeatedl 
taking a large number of 6 values at Magara 
oc eaig The lowest value of Q-'=6/n 
measurable was 5x10~° (This is typ; l 
thiche ypical of most 
Temperature stability was found to be +0.2 
600 K and + 2 K thereafter. ony 
Further, a complete temperature run was made on 
copper of commercial grade purity. The sample was a 
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Fig. 5—Damping versus temperature for polycrystal- 
line commercial grade copper [Frequency of 
oscillation 1 Hz] 


wire of 0.7 mm diameter x50 mm long and was 
annealed in situ at 400°C for 2 hr. The results are shown 
in Fig. 5. It can be seen that the temperature at which 
the peak occurs is in the range for the grain boundary 
relaxation in Cu. The high temperature rise is due to 


‘the background contribution’. The actual peak is 


obtained by subtracting the estimated background 
damping from the measured total spectrum. 


6 Conclusion 

We have described a very flexible system for 
measuring the internal friction in materials. It can 
handle a wide range of frequencies and temperatures. 
An experiment in the low frequency, high temperature 
region has demonstrated the capabilities of the system. 
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The exchange repulsion between ions in ionic solids is discussed. 
Analytical form of the orbitals constructed by Clementi suggests 
many exponential terms which constitute the interacting potential. 
An approximate form is suggested and limitations are discussed. 


The Born-Mayer exponential form (Ae~"”) for the 
repulsive energy between closed shells has been used in 
the study of properties of ionic solids'’*. Here p is a 
parameter. Although simple in form and supported by 
the quantum mechanical calculations* °, it is 
empirical in nature®’!®°. It has a limitation for small 
values of r, the distance between the ions. In fact, the 
Born repulsive interaction may consist of several terms 
partly attractive and partly repulsive’®. Unsold! and 
Bruck'? have shown that the energy has a form 
m(r)e ” where x(r) is a polynomial and y is a 
parameter. Similar is the nature of the expression given 
by Heitler and London*. The Born-Mayer form can be 
taken as an approximation to this with z(r) a constant. 

Recently, Gordon and Kim’* have given a priori 
calculation of intermolecular forces between closed 
shell molecules. Their approach predicts remarkably 
well the distance (c) at which the potential crosses zero 
as well as the distance between nuclei at the potential 
minimum, in case of rare gases. The interaction energy 


between ions of the alkali halides is expressed as [Eq. 


(17) of Gordon and Kim'*]: 


a” i riser iar: | anr3p,(r;) par)! Al) 
aa eget A 1 
where "ie = 1Mz,~1) R + F(R,r;,r2) 
_ 245 _ oe 
(Z,+)R+r; +|R—r,| 
2Z 1 


— @,-1Rtr, +1R—"11 


ee 
FR Ii72)=Rop,4[R—71 1% r,<|R—r;| 


et) 2. ea 
=5\7¢7,) 4rr2  4nrR 


i 
for [R—r,|<r2<IR+n11 


a for r2>R+r; 


where R is the distance between the ions a and b, r, and 
r, refer to the distances from them, Z, and Z, their 
atomic numbers. The electronic charge densities (p, 
and /P,), of the ions can be obtained in analytical form 
by use of atomic orbitals constructed by Clementi!* 
which are given as 


Wiaa= Y Xpaa Ciap ... (2) 
P 


where p refers to the pth basis function of symmetry A(A 
corresponds to quantum number /), « indicates the 
subspecies and i refers to the orbital symmetry A, and 


Xapa = R,,(r) Y ,.(9,~) (0.9) 


where 
Rp=[2(m,,)!]- '/7(2€,,)"4"* */? 


x Pp} e~Eap- 


Here, r refers to the distance, n,, is an integer >A + 1 and 
Y,, are the spherical hermonics. The exponent ¢,, has 
been chosen so as to give the best energy using process 
of optimization. The quantities C;,,, nz, and ¢,, are 
numbers and can be obtained from the Tables of 
atomic wavefunction’*. 

Use of Eq. (1) coupled with Eq. (2) leads to a sum of 
large number of exponential terms for the interaction 
energy which can be expressed as 

a Air 


A; + y Bye 47 
i= r j= 


R 
P — Ayn —Ag 
ao +. O,re 


ee) fe) 
where A,, Bj, C,, Og, Ai, Aj, Ax, As are constants; m, Q, P, R 
are the integers up to which the series is to be summed 
up. They are related to the number of basis functions 
used by Clementi'*. 

It may be noted that Eq. (3) contains a term of the 
type e “/r which may be essential to account for 
infinite repulsive energy in case the ions are brought 
sufficiently near. This type of potential was first 
suggested by Hellman'* for metals and later used by 
Varshni and Shukla!® for alkali hydrides. Puri and 
Srivastava!” have used this form for alkali halides. In 
fact, a form of the type 


Ae~"/r+Be-"+Dre-* ... (4) 
may be a better approximation but such a form is 
linked with the difficulty of evaluating the parameters 
A, B, D, «, B and 6. If equilibrium conditions are used 
for the purpose, one is constrained to limit the form to 
one term only, either the first or the second term of Eq. 
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(4), as the third term alone is unrealistic. Since in most 
computations, one is concerned with the properties in 
equilibrium state, the first or the second term of Eq. (4), 
which are respectively screened Coulomb potential 
(SCP) of Hellman'* and Born-Mayer’s (B-M) form, 
may yield good results'”’'* in case of alkali halides. But 
wherever distances lower than the equilibrium distance 
between the ions such as that in the computation of ‘a 
(Table 1) at which the potential crosses zero, are 
involved, SCP form may show better performance as 
the B-M form has a limitation for small values of r!°!? 

Also, in case of alkali and heavy metal halide 
molecules, results of dissociation energy(Table 2) with 
this form (SCP) agree well with the sophisticated 
calculations of Matcha?°, who has used linear 
combination of atomic orbitals wavefunctions for 
alkali halide molecules. Polarization functions with 
suitable angular behaviour to describe the distortions 
accompanying molecule formation have been included 
in the basis set which comprises atomic orbitals. 

Further, the SCP form gives a value of the Schottky 
defect formation energy for KBr greater than that for 
KCl. This observation is in agreement with the 
experimental observation?! (Table 3). The results of 
calculations of Boswarva and Lidiard2? and Rao and 
Rao’* through the Born-Mayer form are in the reverse 
order. 

Goodisman** has suggested the division of the 
potential for the whole range of r into three groups: 


(a)r<a 
(b)r>a 


a 


Table 1—Calculated Values of ¢ (in au) of Ionic Pairs 


System SCP B-M Gordon and Kim'? 
NaCl 2.05 2.01 2.06 
KCl 2.33 2.39 2.27 


ee 


Table 2—Values of Dissociation Energy 


[in keal mol] 
Molecule Matcha?° SCP 
Lic 109.58 104.77 
LiBr 95.97 90.71 
NaF 111.65 107.5 
NaCl 91.58 93.32 
KF 110,39 115.62 
KCl 95.34 98.80 


Table 3—Energy (in eV) of Schottky Defect Formation 
Crystal Exptl scp B-L R-R 
(Ref. 21] (Ref. 22] (Ref. 23] 
KCl 2.22 2.09 2.05 2.11 
KBr 2.53 2.24 1.92 1.87 
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(c) intermediate range, where a takes certain suitable 
values‘of r. Although we believe that such an approach 
will be quite appropriate, there 1s difficulty of joining 
the potentials at the boundary and it may not always 
be easy. ty Se 

Thomas-Fermi-Dirac method, which is an appro- 
ximate version of density functional method, has also 
been applied’° in case of cohesive energy of KBr. The 
method has not given results. 

While using the commonly prevalent forms of 
potential to cesium and lithium halide crystals, a point 
to be noted is that these metals are the two extremes of 
alkali group; as such their electronic configurations are 
also of two extremes. Whereas lithium has the first s- 
closed shell unshielded by p-electrons, cesium has 
many s and p closed shells. Therefore, their behaviours 
are expected to be somewhat different from those of 
other alkali group metals. Phonon spectra of lithium?® 
and thermoelectric power of cesium?’ are typical 
examples. In case of cesium compounds, the situation is 
complicated on account of the closed shell electron 
wavefunctions ranging to distances much beyond the 
equilibrium distance of halide crystals. 

The penetration of wavefunction into the regions of 
halide ions modifies the distribution of electrons giving 
rise to polarization effects. The effect is increased when 
the halide ions have also long ranging wavefunction of 
closed shell electrons such as iodine. In other words, 
the overlap interaction in case of cesium iodide will be 
the longest amongst alkali halide crystals and in sucha 
situation no single potential form can be expected to 
give good results. 
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Fluorescent wavelength shifting by cold emission of phosphors 
exhibiting ‘quantum-counter-action’ (QCA), popularly called 
summation of photons, has been proposed to enhance the 
conversion efficiency of silicon solar cells for AM1 use. The absolute 
quantum efficiency (number of visible photons produced divided by 
the number of IR photons absorbed) and the optical efficiency (light 
transmitted into the solar cell compared to light incident on the 
fluorescent sheet) of terbium-doped lanthanum fluoride and 
thulium-doped calcium tungstate have been measured. It has been 
verified that there is an optimum screen thickness of the phosphors 
which is a function of the incident wavelength. 


The photovoltaic effect was first reported by Becquerel 
in 1839 and only in 1954, an RCA group! claimed a 
conversion efficiency of 6%. In the later years, 
vigorous efforts have been made to use different 
materials like silicon, gallium arsenide, indium 
phosphate, gallium phosphide, aluminium arsenide 
and even some organic substances? to achieve higher 
conversion efficiency. Most of the photovoltaic work 
has been done in the visible region of the solar 
spectrum. Loferski’s formula? for the maximum 
efficiency of the solar cell indicates that most of the 
visible photons should be used up to produce electron- 
hole pairs and that the IR beyond 1.0 um Only heats the 
solar cells. The solar cell efficiency decreases with the 
increasing temperature*, leading to the loss of 
photovoltage to the extent of 2mV/°C. It is with this in 
mind that the author thought of using ‘quantum- 
counter-action’ (QCA) of a few prospective phosphors 
as efficient solar cells. - 

Literature shows that Hovel® has enhanced the 
conversion efficiency by 2% by fluorescent wavelength 
shifting on the solar cells, using commercially available 
plexiglass. However, there is an Overlap between the 
absorption and emission bands*resulting in a little re- 
absorption of the emitted radiation energy, i.e. a sort 
of self-absorption. In order to look for an up- 
conversion phosphor having a reasonable high 
quantum efficiency (number of visible photons 
produced divided by the number of IR photons 
absorbed) and also high optical efficiency (light 
transmitted into the solar cell compared to light 
incident on the fluorescent sheet), absolute values of 
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these parameters were measured using Bruner’s 
method®, with a few experimental modifications. 

Terbium-doped lanthanum fluoride was prepared 
from spectroscopically pure oxide obtained from M/s 
Johnson Mathey Chemical Ltd, UK. Varying 
thicknesses of thin film of this phosphor were 
deposited on microscope cover slides used as 
substrates. Both vacuum evaporation and chemical 
deposition techniques were employed, and the values 
of areal density were computed® and the quantum 
efficiency calculated in absolute terms at various IR 
radiations along with optical efficiency. IR lamp was 
used with interference filters purchased from Oriel 
Corp., USA. Similar observations were made for 
thulium-doped calcium tungstate phosphors (Table 1). 
The excitation spectra of these phosphors were taken 
with various concentrations of thulium. 

A Reeder thermopile made of gold, purchased from 
Charle’s M Reeder and Co., Michigan, USA, and 
blackened with evaporated colloidal gold, was used to 
measure the incoming IR. The details of this 
thermopile are given elsewhere’. A _ specially 


Table 1—Values of Areal Density, Absolute Quantum 
Efficiency and Optical Efficiency 


Sample Areal Absolute Optical 
number density quantum efficiency 
mg/cm? efficiency 
Terbium-doped (0.05 % by weight) lanthanum fluoride 
I 4.33 0.065 0.055 
2 3.21 0.062 0.058 
3 7.05 0.053 0.052 
4 1.54 0.060 0.032* 
5 0.79 0.054 0.042 
6 2.42 0.048 0.038 
7 0.042 0.064 0.052: 
8 0.850 0.058 0.048 
9 1.30 0.074 0.063 
10 0.520 0.083 0.074 
Thulium-doped (0.06 % by weight) CaWO, 
I 3.48 ~ 0.066 0.052 
2 1.62 0.048 0.016t 
3 5.46 0.054 0.042 
4 3.84 0.048 0.044 
5 2.54 0.062 0.054 
6 0.79 0.057 0.043 
7 0.42 0.064 0.062 
8 1.52 0.074 0.068 
9 1.32 0.086 0.074 
10 0.56 0.065 0.060 
“Looks good for QCA on solar cell 


t Prospective candidate for QCA on solar cell 


(“saw oO 


constructed RCA 931A Photomultiplier-type photo- 
meter was used to measure the intensity of the 
luminescent flux coming from the back of the 
phosphor Screen. The inter-calibration was accom- 
plished in such a way that the final results were 
independent of the absolute intensity of either 
mstrument, i.e. the Reeder thermopile and_ the 
photomultiplier-type photometer. Appropriate cor- 


rections for the absorption and scattering of light 
within the sample were made. 


Both the phosphors exhibited QCA to convert IR 
into the visible luminescence as a function of sensitizer 
concentration for constant condition of excitation. 
The energy transfer in such trivalent rare earth ions 
takes place either by one or more transfers because the 
multiplets associated with these are spaced con- 
veniently and a number of these have long lifetimes 
and show strong fluorescence. Concentration 

_ quenching has also been observed. The absolute 
antum efficiency is not very high, but the optical 
sncy is not discouragingly low either. The 
stence of an optimum thickness is an encouraging 
gn for its use on the solar cells. The IR has more than 
6 of the solar energy at AMI. Its conversion into 
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a J _ > - - es * 
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NOTES 


luminescent visible radiation has two-fold advantages 
and is worth trying on the silicon solar cells. In 
addition to their quantum and optical efficiencies, the 
absorption coefficients for the IR, direct visible flux 
and the visible radiations produced as a result of QCA, 
refractive indexes and also any chemical aging effects, 
etc., are to be thoroughly investigated. 

The work was carried out at the Laboratory for 
Atmospheric and Space Physics (LASP), University of 
Colorado, Boulder, Colorado, USA, and the funds 
were provided partly by NASA. The author wishes to 
place on record his thankfulness to the Director and 
the Dean, IAT, Pune, for publishing this work. 
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das the difference 
F,.6H,O, toa single heat pulse. AT is determine 
ian dearer 6 ego 7 between the equilibrium temperatures of the ee 2 


MnSiF,.6H,O0 and NiSiF ¢.6H20 = the beginning and at the end of the heat paar 
the Phase Transition Points equilibrium temperatures are calculated from*: 


if 


M P SINHA* In(T— T,)=In(T" — T,)— Bt (I) 
“ty aged age epee pny Ba a ass icons where 7” is the instantaneous temperature of the — 
4 calorimeter just at the end of the heat pulse, T, 
em Ai temperature of the radiation shield and B~' tl 
dies taken by the system for attaining the equilibr ri 


E Indian Institute of Technology, bt a seae 
Cryogenic sens aoe aie < =f temperature. The final banger ion tempe erat 
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‘ recognized as the value of 7 at time t=," ze 


The specific heats of CoSiF,.6H,O, MnSiF,.6H,O and plot of 7 versus t and log(T—T,) versus ¢ in t th 
NiSiF;.6H,O are measured in the temperature region 100-300 K rating and in the after-rating pk ed be equ 


using a pulse heating technique. In all the three compounds it is temperatures were determined?. The re 


_ found that specific heat has anomalous temperature dependence 
around the structural transition point. The data are analyzed using the experimental data and the 


a the asymptotic power law. 
dependence of Coan page 


© results of a study of temperature dependence of chloride through the transition 
c heat around the structural transition point!~3 _ tal data agreed within 1 % of t 
rate -d cobalt, manganese and nickel Ne Sige Specific heats C, in the t 
300K for CoSiF ,.6H 20 
7 NiSiF 6. —— are oS Me 


transition peut ‘eds been saccades 
little information‘ is available about the 

fi ¢ heat in the critical region where the 
Fone cone the other. 
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NOTES 


Table |—Critical Parameters in the Phase Transition Region for 7 > T, 

Name of compound T., K a Range of agreement* 4 B E 
CoSiF,.6H,O 259.02 0.34 6x10 7->%t>7x 1074 37.15 64.7 pee 
MnSiF..6H,O 230.5 0.41 1x 10-*>1>2x 109 3.75 10.12 11.52 
NiSiF,.6H,O 234.15 0.32 2x10>1>2x 104 2.46 34.4] 43.4 


* 
Range of agreement between the experimental data and calculated C, from Eq. (2) 


two transition points in CoSiF..6H,O at critical 
temperatures 7.(I)~231K and T.(II)~259K, re- 
spectively. The transition temperature T.(II) is well 
within the limit of error with the value determined 
earlier in our laboratory using continuous heating 
technique’’. T,(II) also coincides with the transition 
temperature 7.(h) obtained from magnetic anisotropy 
data during heating cycle. Figs 1(b) and 1(c) indicate 
that in MnSiF,.6H,O and in NiSiF,.6H,O the 
transition temperatures are 230 and 234 K respectively. 
The transition point in MnSiF,.6H,O agrees well 
within the limit of error to the transition temperature 
T{h) obtained from the optical absorption data”. The 
transition in NiSiF..6H,O is similar to the magnetic 
transition observed in Ni(ClO,).6H,O [Refs 12, 13]. 
These confirm the presence of hysteresis effect during 
the structural transitions in CoSiF,.6H,O and 
MnSiF,.6H,O. A comparison of the specific heat data 
for the compounds reveals that the shape of the curve 
around the transition point in MnSiF,.6H,O, 
NiSiF,.6H,O and around 7.(II) in CoSiF,.6H,O has 
characteristic similarities. In the region JT<T,, the 
specific heat rises sharply with temperature whereas in 
the region 7>T,, the specific heat decreases at a 
comparatively slower rate with temperature. The 
shape of the curve in the critical region resembles a 
mirror image of J-shaped C, versus T curve obtained 
near the critical temperature in magnetic 
transitions'*"'°. 

Recent theoretical advances in the description of 
critical behaviour within the critical region suggest 
that the specific heat of a system is expected to obey the 
asymptotic form'®, 


C,=At *+B+Et 


} 


. (2) 


| where t =(T—T,)/T., and A, #, Band E are constants to 
be determined from experimental data. The first 


asymptotic term in Eq. (2) corresponds to a singularity 
at T=T, and other terms represent the background 
contribution and are assumed to be linear in 
 temperature’’. The critical region parameters 
determined from the experimental data are shown in 
1. It is worthwhile to note that the parameter 
for T>T, only. At T<f,, the variation of 
with temperature is very sharp and it has 

sen possible to obtain more than one point in this 


(a) CoSi Fg6H,0 


(b)Mn S: Fe 6H, 0 

(c)Ni Si FS 6H,O 
ee Experimental poms ‘ 
‘7 Calculated points 

ey from £q.(2) 


© Experimental 
points 
\, ~~” Calculated po 
trom €q(2) 


Colcal g-mot” K') 


0-001 0-01 0.000! 
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Fig. 2—Specific heat of (a)CoSiF..6H,0,(b) MnSiF..6H,O and (c) 
NiSiF..6H,O as a function of reduced temperature f in the vicinity of 
the phase transition for T> T, 


region. Table 1 also shows the calculated critical 
temperature and the ranges of agreement between the 
experimental data and the specific heat C, calculated 
from Eq. (2). Fig. 2 shows the variation of C, with 
reduced temperature f plotted in a log scale in the 
transition region. Experimental points when joined 
together indicate a linear variation while calculated 
curve shows deviation from linearity. This deviation 
probably arises from the assumption of linear 
temperature dependence of the background contri- 
bution in Eq. (2). The experimental plot in Fig. 2 does 
not show any rounding off near the transition 
temperature and for small positive values of the critical 
exponent «, the logarithmic singularity of specific heat 
is obtained for each sample. The specific heat in the 
vicinity of first transition 7,(1) in CoSiF ..6H,O has not 
been analyzed using Eq. (2) because the experimental 
data do not give meaningful variation of C, versus log t 
and cause a large percentage error in the calculation of 
critical parameters. 

In conclusion we give an evidence for logarithmic 
singularity in the specific heat of paramagnetic 
compounds around a structural transition point. A 
similar observation has been made by Sinha et al.* in 
case of the other paramagnetic compounds around 
their structural transition points. It is presumed that 
the present results support that structural transition 
has a second order origin. 
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A fresh ‘study of the harmonic potential constants, compliance 
_ constants and vibrational mean amplitudes has been attempted in 
some pyramidal X YZ, cases (C,) using the kinetic constants method. 


The kinetic constants’ * relating to molecules have 
been studied at length and these provide interesting 
results concerning molecular dynamics. The kinetic 
constants have been advantageously utilized in solving 
the secular equation, governing potential constants. 
The molecules of pyramidal XYZ,-type belongs to C, 
point group having six fundamental frequencies under 
two species (4A’ + 2A”). The present investigation has 
n undertaken to obtain the reasonable set of 
a, constants for some pyramidal XYZ,-type 
eri the most general quadratic valence 
1. In a addition to potential constants, other 
constants such as me a constants 
rac mean ‘amp itudes’ of bonded « 


quadratic potential energy function has _ been 
considered and the F matrix elements have been 
obtained as follows: 


A’ type: 


So 4fop 
d? fy 


V2 Soa 
J 2d fap 
Sat Saa 


\/2Da fos 
of: 2Dd* fu 
vba fae +5 + fix) 

Dd ( a+ 


A” type: 
( ate va Ghee fk ) 


The elements of kinetic energy matrix’ are used to 
calculate the kinetic constants of these molecules. 

The method of kinetic constants seeks to relate the 
off-diagonal elements to the diagonal elements of the 
F-matrix through the relation: 

Fj, ki . A eek Brel oy 

cL, (i< ji i, j=l, 2,3, 4) 

The compliance constants are also per by the 
Decius method’. 

Using Cyvin’s equation 2 = LAL, the symmetrized 
mean square amplitudes and hence the valence mean 
square amplitude ae for both the cee i: 
the hs foliea distane ated at 29 
using the present set of 
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Table 3—Values of Compliance Constants (A/m dynes) 


nN, Np 


Molecule Np Ng Nga —Npg 
AsHD, 0.3271 0.5346 0.0381 0.0003 0.7259 0.2526 
AsDH, 0.3270 0.5077 0.1280 0.0003 0.8633 0.7545 
AsHT, 0.3272 0.5755 0.2739 0.0001 0.8285 0.1246 
AsTH, 0.2893 0.7451 0.3851 0.0001 0.5626 0.4305 
AsDT, 0.3101 0.4577 0.1193 0.0003 0.7884 0.2979 
AsTD, 0.2929 0.3852 0.0303 0.0005 0.7192 0.4959 


Nay — Nyy — Np, Npp — Nyy Ny — Nag 
0.1162 0.0003 0.0036 0.0002 0.0021 0.0001 0.0077 
0.2040 0.0043 0.0037 0.0002 0.0023 0.0015 0.0101 
0.3148 0.0006 0.0022 0.0001 0.0105 0.0052 0.0048 
0.3562 0.0038 0.0020 90.0001 0.0047 0.0025 0.0057 
0.1512 0.0013 0.0035 0.0002 0.0082 0.0023 0.0151 
0.111% 0.0006 0.0038 0.0003 0.0022 0.0021 0.0151 


Table 4— Values of Mean Amplitudes of Vibration (107 2A) 


298. 16K 

Molecule lp ly I, l, 

(X—Y) (X—Z) (¥..;2) (2:22) 
AsHD, 8.8533 7.2927 15.1453 12.7990 
AsDH, 7.9362 8.8430 14.7773 14.7781 
AsHT, 8.7714 6.5792 14.2793 11.9297 
AsTH, 6.6569 8.3665 13.8342 15.5342 
AsDT, 7.8348 6.8434 12.7785 12.6163 
AsTD, 6.6171 7.7660 12.7546 12.9721 


the major stretching kinetic constants kp of As—H is 
half that of As—D and one third that of As—T. The 
stretch-stretch interaction kinetic constants kp, and k, 
practically remain the same for the isotopic pairs. The 
interaction kinetic constants, viz. kpp and ki, are 
uniquely negative in all the cases. 


As the mass of the Z atom increases, the major 
stretching force constant fp increases. This trend can 
be noticed in ASHD,, AsHT,; AsDH,, AsDT, and 
AsTH,, AsTD). The stretch-stretch interaction force 
constant /,4, angle-angle interaction force constant Say 
and the bond-angle interaction force constants Jog and 
Jax ASSume negative sign in all these cases. 


The compliance constants for the molecules under 
consideration are listed in Table 3. These constants are 
invarient to the choice of coordinates defining the force 
field and they may be used as a measure of the bond 
strengths and interactions instead of force constants as 
pointed out by Decius® and Jones!2. Further, it may be 
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noticed that the compliance constants exhibit trends 
opposite to that of the force constants. 

The vibrational mean amplitudes for both the 
bonded and non-bonded distances at 298.16K 
reported in Table 4 are in the expected range. It may be 
noted that /p remains same in isotopically substituted 
molecules, viz. ASHD,, ASHT,; AsDH,, AsDT, and 
AsTH,, AsTD). It may be added that a decreasing 
trend is noticed in /, when the mass of the Z atom 
increases as per the order cited. 

The authors are thankful to Prof. S Sathikh, 
Director, Madras Institute of Technology, for his kind 
and constant encouragement and facilities given to 
carry Out this research work. 
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Super position of Potentials for Calculation of x, and w,.x,—The potential functions 
Diatomic Molecules (1) and (2) were differentiated up to fourth degree which 
S V BIRAJDAR is essential for the calculation of %, and w,x,. In both 


cases the first derivative at r=r, iS zero, which is in 
lege of Science, bad 001 : , ae 
SB Coleg of Rerenpeond 3 accordance with the required criterion for extremum. 


P L SARDESAI & S H BEHERE* The successive values of differentials at r=r, are as 


follows: 
Department of Physics, Marathwada University, Aurangabad ; : 
? 431 004 : 2} Potential (1) Potential (2) 


Recejved 16 August 1982 ul (r.) a 9 7 p2 Fe3b? 

The potential functions for diatomic molecules suggested by lyer pe fe aS... Seg 
and Sharma [Indian J Pure & Appl Phys, 17 (1979) 615] are U(r.) ae =b or i 8b 
reinvestigated. A new potential function is also suggested. The «, and 
x, values of a few diatomic molecules are calculated using the u(r.) D. 227 44 D. 17b4 
newly suggested potential function and these values are compared . 2-15 ae 
with the experimental values. The new function gives better Sea Sa 
agreement as compared “ieaaaen ones ve a bylyer Following Varshini*, 
and Sharma. 


? Xr, 6 B2 
. The importance of potential functions for the diatomic %e= | 3 | o, 
molecule is well known. Various molecular constants . 


a. i ti ee) -” 


4 _ like dissociation energy and anharmonicity constants where X =U") esr . .. (3) = 
a ean be calculated using the potential functions. = ™ ve = J 


 & ne ions ethics ane superimpositions of potential 
on ela tution se The sees ot ss 
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Table 1—Values of x, for Different Diatomic Molecules 

Molecule a, (Exptl ) Potential (1) Potential 2) . Present Potential (5) 

mr ae ES "Error "a oe Error a, Baas 

(Calc) oe (Calc.) Me (Calc.) ei 

Br, 0.0003187 0.0005629 + 76.0 0.0004470 + 40.26 oe 60 i ve 
CH 0,534 0.7496 +40.4 0.5493 +286 04 So ee 
Cl 0.00149 0.002429 +735 0.001912 +36.59 0.00176 ey 
CO 0.01750 0.02271 +297 0.017146 ..—204.. 0015675 . = 08s 
CdH 0.2180 0.4730 +116.97 0.3708 +70.12 0.3438 te 
H, 3.0620 2.0824 — 31.99 0.8799 —71.26 0.5623 —8l. 
HBr 0.2333 0.3504 +50.23 0.26398 +13.16 0.2411 + 3.34 
HCl 0.3072 0.4364 +4207 0.3248 + 5.7 0.2954 — 3.84 
HF 0.7980 0.9222 +15.56 0.6596 —17.33 0.5899 — 26.08 
HI 0.1688 0.2659 +57.53 0.2030 +20.24 0.1863 + 10.39 
HgH* 0.2060 0.2678 +30.02 0.2008 —2.5 0.1931 — 6.26 
I, 0.0001138 0.0001994 +75.20  0.0001592 + 39.88 0.0001486 + 30.55 
ICI 0.0005354 0.0008 746 +63.34 0.0006928 +29.39 0.0006447 + 20.4 
Li 0.007040 0.01435 +103.0 0.01063 + 50.95 0.009646 + 37.02 
K, 0.0001650 0.000563 1 +241.0 0.0004327 + 162.25 0.0003983 + 141.38 
N, 0.01732 0.02685 +55.02 0.02060 +18.99 0.01896 +9.4 
NO 0.01710 0.02552 +49.21 0.01973 +15.36 0.01819 + 6.4 
Na, 0.0008736 0.002034 +133.0 0.001539 +76.17 0.001408 +61.19 
OH 0.7242 1.0027 +38.45 0.7368 + 1.74 0.6667 -—79 
O, 0.01593 0.02345 + 47.5 0.01816 +1421 0.01676 +5.4 
P, 0.00149 0.002351 +57.82 0.001836 + 23.2 0.001699 + 14.00 
SO 0.005736 0.008115 +41.47 0.006289 + 9.6 0.005807 + 1.20 
ZnH 0.2500 0.5638 +125 0.4406 +76.25 0.4081 + 63.20 
AuAl 0.000668 0.0007488 +12.09 0.,0005713 —14.52 0.000524 —21.55 
GeSe 0.0002890 0.0004 146 +43.47 0.0003217 +11.33 0.0002972 +2.8 
He} 0.2240 0.3242 +44.76 0.2386 +6.5 0.2159 — 3.62 
Si, 0.001350 0.001896 +40.43 0.001463 + 8.3 0.001349 — 00.09 


Sc a ee 


potential functions can be formed as 
D, 2 
Ur)=5* {[1 —exp(—ap)] 


—[1+bp] exp(—bp)+1} Fae 
where a and b are Morse’s and Rydberg’s parameters 


and bP a Ke which means 


e 


; K 
+ e 
respectively, and a =sp 


e 


— b — 
that may or b=a,/2. 


The above suggested potential function can be written 
either in Morse’s or in Rydberg’s form. 


In Morse’s form 
D, 4 
U(r)=—* {Ll —exp(—ap)] 


—[1+a,/2 plexp(—./2ap)+1} (6) 
In Rydberg’s form 


wt fee (- 


~ (1 +bp] exp (—bp)+1} - (7) 


Potential function (5) satisfies all the required 
criteria. Atr=r,, in the potential functio ns (5), (6), (7), 
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U(r.) =0 as the last term is a whereas in the potential 
functions (1) and (2), U(r.)= —D, as the last term is 


D. : 
THe This has no effect on the evaluation of a“, and 
@,.X- which need the successive differentiations and the 
last term being only an adjustable constant its 
derivative is zero. The values of successive differentials 
at r=r, are: 


Morse’s form Rydberg’s form 


D 


U" (r.) —.4a, 5 2? 
D.| 3 D.| 3 
Up, cas 3 bod 2| 3 oe 3 
(r,) 5 E 4a ; = +2] 
D D. 13 
U'(r.) —.26a* — = 
2 aa 


By calculating the values of X and Y and using 
relations (3) and (4), the values of a, and @,X, are 
evaluated. Although Eqs. (6) and (7) look different, they 
are basically the same; naturally the relations for «, and 
®X¢ derived either by Eq. (6) or Eq. (7) are the same. 
Thus the potential function (5) gives ae 


%=(0.9714 Az 1) 9Be 


ew Ts 


NOTES 
j Table 2—Values of w,x, for Different Diatomic Molecules 
F Molecule wx, Potential ay Potential (2) Potential (5) 
(Expu ) OX, Error @X, Error WX, Error 
(Cale.) % (Calc.) _ (Cale.) ¥ 
Br, 1.0774 2.398 +122.0 1.69 +570 1.57 + 46.0 
CH 63.0 101,37 +60.90- 71.76 +1390 66.61 +5,70 
Cl, 2.67 5.39 +102.00 3.82 +43.00 3.54 + 32.0 
co 13.28 18.92 +4200 13.40 + 1,00 12.43 ~ 6.00 
CdH 46.30 106.00 +128.00 75.09 +62.00 69.63 + 50.00 
H, 121.33 182.00 +49.00 130.15 +7,20 120.80 -0.2 
HBr 45.22 80.73 +78.00 57.17 +26.00 53.00 +17.00 
HCl 52.82 87.47 +65.00 61.90 +17.28 56.80 + 8,00 
HF 89.88 120.06 +33.58 85.02 — 5.40 78.83 ~ 12.83 
HI 39.64 75.64 +90.0 53.56 +35.0 49.67 +25.0 
HgH* 40.90 59.55 +45.60 42.16 + 3.07 39.10 —4.38 
: a 0.614 1.33 +117.0 0.0450 +52.0 0.876 +43.0 
Sart “Gore lao UC} 1.501 3.06 +1040 216 + 44.0 2.01 + 34.0 
ger Li, 2.610. 5.25 +101.0 3.717 + 42.0 3.44 + 32.0 
2 a 0.2829 0.7354 +159.9 0.5205 +84.0 0.4832 +70.8 
N, 14.32 25.34 +77.0 17.94 | +250 ~ 16.64 + 16.0 
NO 14.08 24.72 eroOs”--* ESI +24.4 (G29 +150 
Na, —_—0.7254 1.56 A2kS.0s% 1.108 +528 1028 +420 
—— bei Os: (a. 84,88 136.36 = +. 60.6 96 56 os 413.7 89.54 +5.4 
: Ogpsices WOR, a Gh cpt 806+ 1530. +277. 1449 * +184 
a Peete eee oa. EGER» 2.64 7 2 26S > Fas —15.0 
ae SEN; 0: ., SRB cy, 18 ORO y 2. TL PEs = SRR FS © 2-720 + 28.0 
oe Zn SSF 12281 4122.0 86.97 +57.0 80.6 +46.0 Ts <3 
a AvAl —s-1.1630 1.4891 = +:28.03 1.0541 - 9.6 0.9784. = 15.86 7 ats 
GeSe _—_ 11.36 1.5052 +10.68 1.066 —21.62 0.9896 —37.43 | , 
35.30 52.70 +493 37.31 + 5.69 34.61 — 1.95 
2.02 3.34 —+65.40 2.36 +1708 219 +8.68 


een Table 3—Values of A for Different Diatomic Molecules = ; aig 
 . - ; , : J 5 ¢ Se Nr OLD oa ; Dig 6 iF : ‘ V \ : . - 
- _ Molecule State iy’ i J ts i A > re, Az sath re “a =P, 35 Bit yr Sa tees DR, rae i a ~ : ae 
Sa (Varshini) (Present a tee a ag a al ae 


~ es 
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ee 7.655 AW .. 9) 
La r e 

The values of x, and @,x, using relations (8) and (9), i.e. 

potential (5) are calculated and compared with the 

experimental values (Tables | and 2). 

Discussion—It was needed to revise the values of A 
for the molecules under consideration as very recent 
and authentic data> of molecular constants compared 
with those of Varshini* had become available. The 
values of A obtained by the present calculations are 
close to those of Varshini (Table 3) except for SO. 
Instead of N,(ii), CO(ii), Cd(iii) and NO(ii) reported by 
Varshini*, we have included a few other molecules, viz. 
GeSe, AuAl, Si, and He}. Instead of HgH, HgH* 
molecule is selected as the wx, and &, are to be derived 


eS: from interpolation of AG, and B, values. These 


alterations change the average % errors in %, and ,x, 
__ using Morse and Rydberg potentials, which are not 
_ explicitly reported in this note for the sake of brevity. 
The average % errors in , for Morse and Rydberg 
functions using these randomly selected molecules are 
25.98 and +22.39 respectively (Varshini’s Table 
a Ay. . The average % errors in %, using potentials (1) and 
= 2 ane a + 66.73 and. +32.33 respectively, whereas it is 


+ 26.54 only when potential (5) is used (Table Dw 
is in close agreement with the values derived i 
Morse and Rydberg potentials. 

In the case oj w,x, (Table 2), the average % erTO1 
Morse and Rydberg potential functions using revi 
are +26.82 and +19.92 respectively; whereas t 
average percentage errors using potentials (1) a and 
and (5) are + 81.12, + 30.48 and +24.48 respe ective 


Conclusions 
It can be concluded that the superposition of | 


and Rydberg potential functions suggested ~ 
and Sharma has not given very rem 
improvement in the calculations of , < 


Sharma. 
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The screened p-p reaction rates to form *H have been studied 
considering the effect of screening on nuclear reaction rates in stellar 
interiors. The scrreening effect appreciably augments the reaction 
rates at the lower end of hydrogen burning temperature. As the rates 
are greater, mean life of the star in this phase of evolution will be 
shorter. The screened rates give us a new set of data to reinvestigate 
‘ the luminosity, mean life and relative abundance of 7H. 


In the main sequence stars, the thermal pressure due to 
the hydrogen burning reaction establishes a balance 
with the gravitational pressure. A good study on 
various aspects of the reaction such as Coulomb 
: gh penetration, etc was made by Burbidge et al.’ 
_ At a temperature of the order 10°-10’K, hydrogen 
burning becomes dominant. The transmutation of 
hydroge n to deuteron is given by: | 

TH+ BY ty sgl} 

eutrino quickly escapes — Siete the star and 
‘on combines s with an electron to give gamma 


‘ 2 
“gy - ie 
= ts. a 


the pp reaction and 


ct ih neh i | 


plasma, and z, ¢ and k are respectively atomic number, 
electronic charge and Boltzmann’s constant. 


We may write 


x 
n.=pN, - for hydrogen plasma ...(4) 


where p is the density, N , the Avogadro number, X the 
fractional abundance, A = 1 for hydrogen ions. 


Taking pX = 10?gcm > 
n, = 6.025 x 107° 


, We get 
.. (4a) 


We find that T < 1. So the reaction is in weak screening 
regime. 

For nuclei to react, it is necessary that they should be 
able to approach each other closely. If the nuclei are to 
approach a distance r, the kinetic energy for their 


Z122e° 


relative motion must be . The mean thermal 


energy Is = and the energy of relative motion is 3kT 


when they approach (Amburtsumayan"). Hence | 


‘pe (5) 


21230" 


E = 3kT = ———_ r 
r ; 


where r= 107 ve cm, the upper limit at which the nuclei * 


iabeaiekate be 10 em 
partie ee co 1 poss 
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Putting ¥Y =exp(3''*I-*’) and taking logarithms on 
both sides, we find 


0,091 
(8) 
log Y= 0.s(“r a) 


=CT, >! ,««(9) 
where C is a constant. 

Reaction and energy generation rates—The equation 
for unscreened reaction rate (R) is obtained from the 
work of Burbidge er al.’ as obtained in the pape by 
Cohen?. 

ie Z,22M,M, se 2,7-t 
R=7.20 x 10 in| gee 7M, | Sot“e 


= 1568.187 x 10° te" ' .. (10) 


where 


1 =42.48 (Hine To 


M,+M, 
= = 42484751)! (10a) 
ae 10°77 keV barn (Cohen®), ny pelos 


x 1075 for p-p reaction. 


: ‘The logarithm of unscreened rate is 
; eee is 043s ttt) 
|The Sanath of screened rate (R) is 


--(12) 


Table 1—Reaction and Energy Gententeee Rates, as 
Function of Temperature = 


kT T, = 7/10° logR  logR* loge 
2.0658 
4.2782 
5.4522 
6.2068 
6.7340 


7.14459 


1.4147 
4.048 
5.327 
6.1247 
6.6758 
7.1002 


0.1 1.16 
0.2 2.32 
0.3 3.48 
0.4 4.64 
0.5 5.80 
0.6 6.96 


— 2.3095 
—1,5119 
— 0.9608 
— 0.5364 


lower end of temperature is appreciably high. Atl 
temperature, the screened and unscreened ra 
almost same. We can observe from Eqs (2), (3) 
that c*, 


Med 
A 


Sy = constant 


effect at higher temperature, the density she 
At constant density condition, the 
enhancement factor is given Dy. 


=constant. 107°°* 


The screened rates diffe. 
Burbidge et al.’ at the low 


temperature. The resu 
_ find rnean life, luminosity an a 


eee 


sequence star. The 
alculated _ throug b 


Indian Journal of Pure & Applied Physics 
Vol 21, May 1983, p 319 


B System of SnBr Molecule 
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The spectrum of tin monobromide in the range 3250-2950 A has 
been excited in high frequency discharge and recorded in the plane 
grating spectrograph at a dispersion of 3.6 A/mm. The vibrational 
analysis has been carried oyt and the following constants have been 
evaluated: 

v.= 33064.6 cm! 
B?Z: wi =298.12 cm !; wi xi=2.56 cm! 
XM, 3: = 247.00 ont wx) =0.77 cm! 


The spectrum of tin monobromide is known to be 
consisting of four groups of bands. Sharma and 
Venkateswarlu! excited a group of bands in the range 
6400-5100 A in an uncondensed transformer discharge 
and assigned an electronic transition A? £* > X 711, 5. 
_ Jevons and Bashford? excited the spectrum in heavy 


se - = 


Js, viz. 4255-3709 A, 3428-3021 A, around 3600A 

30-2850 A. Oldershaw and Robinson? reported 
up of bands below 2500A in absorption. Jevons 
fc led seven bands in the region 


discharge source and reported four groups of 


a “My 2: 


Table |—Band Head Data of B72 X 711, , System of SnBr 
Rel. 


Wavenumber Assignment Av= Wavenumber 
intensity in vacuum, rie Vabae Veal reported by 
cm! cm”! Ref. 2 
cm”! 

l 33953.0 3.0 —0.3 — 

2 33670.4 2.0 —0.2 — 

2 333828 . 1.0 +0.1 a 

2 33136.9 1,1 —0.3 — 

6 33089.8 0.0 0.0 33091.0 

5 33893.0 1.2 —0.3 — 
10 32845.1 0.1 +0.8 32846.5 

8 32600.6 0,2 +0.2 32601.9 

3 32409.6 1.4 —0.5 — 

6 32357.9 0.3 —0.1 32361.2 

3 32170.6 1.5 0.0 -— 

3 32117.4 0.4 +0.2 32121.8 

4 31878.0 0,5 +0.1 31880.0 

3 31640.2 0.6 +0.1 31646.0(?) 
2 31403.4 0.7 05° 

2 31169.5 0.8 +0.3 — 

3 30936.3 0.9 +0.2 _— 
Bashford?, eleven additional bands have been 


observed. The bands have been analyzed by taking the 


(0.0) band at 33089.8 cm~!' and the following - <a a8 
' ie. 


ae we 
Pe a 


vibrational constants have been cvalianess 


M3 3064,6.cm-=) sas 
BE: w= 298.12 em~!s « re ss 
\ oe a 


«9 = 247.00 cm™ pigs aie 


%, 5% os tp 


os £2 ty 1’ 
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Effect of Dielectric Constant on 
EDTA-Metal Chelates—An 
Ultrasonic Study 
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Employing a single-crystal, variable-path interferometer, five 
EDTA-metal chelates, viz. EDTA + CrCl,, EDTA +ZnSO,, EDTA 
+CoCl,, EDTA + SrCl,, and EDTA + KNO, have been studied at 
30°C and pH=3.5 using 90%-10% water-dioxane mixture as 
solvent. The condition for a metal ion to form a chelate with EDTA 
under experimental conditions has been deduced in terms of radius 
of the metallic cation and valency of the cation. The experimental 
results have been discussed on the basis of Bjerrum’s theory of ion 
association. The results on these chelates in water-dioxane mixture 
have been compared with those obtained when water was used as 
solvent. The present studies have been found to be in agreement with 
those obtained in the case of six other EDTA-metal chelates studied 
earlier by the author. Lower dielectric constant of the medium has 
been found to be favourable for the extent and scope of chelation. 


A chelate is a ring structured stable complex ion. 
EDTA is well known for its capacity to form chelates 
with metals. Ultrasonic studies of EDTA-metal 
chelates, though limited, have proved quite interesting 
in view of their wide scope. Various workers! ~* work- 
ed on these lines but their studies have little scope for 
systematization. The author attempted a systematic 
study on these EDTA-metal chelates, especially in 
regard to the influence of various factors like valency 
(z), radius of metallic cation (r), temperature (7), 
dielectric constant (e), etc., on the process of ion- 
association in EDTA-metal chelation. The results 
pertaining to some metals have already been 
published°~*. The present note reports the results 
obtained with some other metallic salts using water- 
dioxane mixture as solvent to Study the influence of low 
dielectric constant on chelation. 


EXPERIMENTAL DETAILS 
Mixtures studied—Using disodium salt of EDTA, 
five mixtures, namely, (1) EDTA +CrCl;, (2) EDTA 
+ ZnCl,,(3) EDTA +CoCl,, (4) EDTA + SrCl, and (5) 
EDTA+KNO, have been Studied. It has been 
established by earlier studies’ that EDTA cannot form 
a chelate with KNO, using water as solvent. In the 
Present investigation, EDTA+KNO, mixture has 
also been included and studied with a view to knowing 
whether the present solvent of less dielectric constant 
can give any scope for EDTA to form a chelate with 
KNO, also. 
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Preparation of solutions—Using double distilled 
water and distilled AnalaR sample of dioxane, a stock 
solution of 90% -10°%, water-dioxane mixture has been 
prepared. This mixture has been used as the solvent in 
the preparation of EDTA-metallic salt mixtures in 
different molar proportions adopting Job’s method? of 
continuous variation. Using the water-dioxane 
mixture as solvent, equimolar stock solutions of 
individual EDTA and metallic salts of 0.1 M each have 
been prepared. The total molarity of the mixture has 
been kept constant at 0.1 M in all cases of metallic salts. 


Temperature and pH—AIl the mixtures have been 
studied at a constant temperature of 30+0.1°C and 
pH =3.5 using a thermostat of NBE type and a pH 
meter of CL 44 type (Toshniwal Co.). The pH has been 
maintained constant at 3.5 for the reasons explained in 
an earlier paper?. 

Velocity determination—A single-crystal, variable- 
path interferometer working at 2 Mc/s (model made by 
NPL, New Delhi) capable of yielding velocity 
measurements to an accuracy of 0.013 %% has been used’ 
for the determination of ultrasonic velocity by the 
usual well-established interferometric method. 


Density determination—A bicapillary pyknometer — 
and an accurate analytical balance capable of 
measuring the mass to an accuracy of 0.0001 2g (0.01%) 
have been used for density determination. 


Calculation of compressibility—Using the values of 
ultrasonic velocity (v) and density (d), the adiabatic 
compressibility (8) in each mixture and at different 
molar proportions has been calculated using the | 
formula B = 1/v2d - a 


RESULTS AND DISCUSSION . : 
It has been observed that the coloured metallic ions” 7 
formed deeply coloured chelates and colourless 
metallic ions colourless chelates at 1:1 molar. ; 
between EDTA and metallic salt. The variatio ew 
adiabatic compressibility with molar fraction of the 
metallic salt in each mixture has been illustrated in 
Fig. 1. Whereas in four mixtures, B deviated from 
linearity having a maximum deviation at 0.5 molar. 
fraction, ie. at 1:1 molar ratio between EDT, 

metallic salt (Fig. 1), in the case of fifth mis 
EDTA + KNOs, it varied li early® with mola 
(6) at cinaich tla alee 
mixture have been noted from the « 


o 


- er 
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Table 1—Values of z/r, r/z and Maximum Extent of 
Chelation in EDTA-Metal Chelates 
Metallic Radius Valency =r += Maximum extent of 
salt of metal of metal A chelation in terms of 
ion (r) A ion compressibility 
(Ref. 10) 10° '* em?/dyn 


In water In water- 
dioxane 
mixture 

CrCl, 0.65 
ZnSO, 0.74 
CoCl, 0.8? 
SrCl, i 
KNO, 1.33 


~ 


—iIyviv iv 


chelation in terms of compressibility are given for 
different mixtures in Table 1. 

Extent of chelation—The effects of r and z on the 
scope of association in general and EDTA-metal 
chelation in particular have been thoroughly discussed » 
in the earlier papers*-°. To compare the effects of r and 
= on the extent of chelation in the present studies with 
those obtained when water was used as solvent, the 
variation of the maximum extent of chelation in terms 
alors pveseibiliy wi ith picand re of the pester 

n 
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IN WATER 

IN WATER DIOXANE 
PREVIOUS SALTS 
PRESENT SALTS 


MAXIMUM EXTENT OF CHELATION IN TERMS OF 
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Fig.4-—Variation of maximum extent of chelation in terms of 
compressibility with r/z of metal ions in water-dioxane mixture and 
water 


Scope of chelation—The variation of maximum 
extent of chelation with z/r of metal ion is the same as 
that obtained with the salts studied previously® (Fig. 3). 
This curve, when extrapolated downwards, cuts the =/r 
axis at 0.83. This is the minimum value of =/r for any 
metal ion to form a chelate with EDTA. Similarly, the 


straight line plot in Fig. 4, representing the variation of | 


maximum extent of chelation with r/z, coincided with 
that obtained with the previous salts (Fig. 4). This 
straight line, when extrapolated downwards, cuts the 
r/z axis at 1.2 A. This is the maximum value of r/z for 
any metal ion to form a chelate with EDTA. As the 
extrapolation of a curve in general involves 
inaccuracies, the optimum value of r/z has been used 
for the discussion of the results from now Onwards. 


The optimum value of r/z in the case of water as 
solvent, as reported in the earlier papers*® is 1.08 A, 
whereas it is 1.2A in the present studies with water- 
dioxane mixture as solvent. These optimum values of 
r/z reveal that (1) metal ions having r/z value greater 
than 1.08 A have no chance for chelation with EDTA if 
water is used as a solvent: (2) if the r/z values of the same 
metal ions are greater than 1.08 A, they can have the 
chance for chelation with EDTA if water-dioxane 
mixture is used as solvent. In the case of EDTA 
+KNO, mixture, the value of r/z for metal ion, viz. 
1.33A, is greater than 1.08 as also 1.2 A. That is why no 
chelation has been observed in either of the dielectric 
media. The cause for the higher optimum value of r/z 
and the resulting wider scope of chelation in the 
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present studies may be attributed to the lower 
dielectric constant of the medium. 7 

It can therefore be concluded that the two possible 
aspects with reference to which the dielectric constant 
can influence chelation unequivocally point out to the 
fact that the chelation increases as the dielectric 
constant of the medium decreases. 

The above conclusion or deduction from the 
experimental observations can also be explained 
theoretically on the basis of Bjerrum’s theory of ion 
association. Bjerrum’s expression'' for the minimum 
distance of approach (q) for the ions to be called as 
associated, namely 


(where the letters have their usual significance) tells 
that when -,,=,,k and 7 are constants, lower dielectric 
constant (€) increases the value of g which in its turn 
gives more scope for ions to associate. It means that the 
scope of association or chelation is inversely related 
with dielectric constant of the medium. The dielectric 
constants’? of 90°.-10%, water-dioxane and that of 
water at 30 C are 76.73 and 67.98 respectively. 

The optimum values of r/= which are measures of 
scopes of chelation when water and water-dioxane are 
used as solvents are 1.08 A and 1.20A respectively. The 
product of the optimum value of r/z in water and 
dielectric constant of water is 1.08 x 76.73 =82.87. In 
the case of water-dioxane mixture, the product is 1.2 
x 67.98 =81.58. The product in the two cases is 
approximately the same. So, these experimental results 
are in accordance with the theoretical deductions 
based on Bjerrum’s theory of ion association, both 
qualitatively and quantitatively. 

The author expresses his thanks to Prof. Bh. Krishna 
Murthy for his valuable suggestions and to the 
authorities of the University Grants Commission, New 
Delhi, for financial help. 
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